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To study the effect of a strong magnetic field B on the sound velocity Vs of plane waves propagating 
in a strongly magnetized quark-gluon plasma (QGP), a chiral magnetohydrodynamical (MHD) de- 
scription of a perfect (non-dissipative) QGP exhibiting dynamical chiral symmetry breaking (D^SB) 
I is developed using the effective action of the Nambu-Jona-Lasinio (NJL) model of QCD at finite 

. temperature, finite baryon chemical potential and in the presence of a strong magnetic field. Here, 

CNJ ' the D^SB arises due to the phenomenon of magnetic catalysis. Apart from an interesting frequency 

dependence, for plane waves propagating in the transverse or longitudinal direction with respect to 
' the B field, the sound velocity is anisotropic and depends on the angle between the corresponding 

, wave vectors and the direction of the B field. Moreover, for plane waves propagating in the trans- 

' verse (longitudinal) direction to the external B field, the sound velocity has a maximum (minimum) 

\ at T < Tc, reaches a local minimum (maximum) at T ~ Tc and remains constant at T > Tc. Here, 

• Tc is the critical temperature of the chiral phase transition. Thus, the constant value ~ 1.5cs at 

' T >Tc turns out to be a lower (upper) bound for waves propagating in the transverse (longitudinal) 

. direction with respect to the external B field. Here, Cs — l/VS is the sound velocity in an ideal gas. 

' PACS numbers: ll.lO.Wx, 12.38.Mh, ll.30.Rd 

T— I ■ I. INTRODUCTION 

> ■ 

I Hydrodynamics is a framework to describe the time evolution of a system under local thermal equilibrium (LTE). 
■ Its relativistic version is used extensively in the past years to describe the phenomena observed in the heavy 
ion experiments at BNL-RHIC, where excited nuclear matter has been already created in the Au-Au collisions 
lO ' with the collision energy of ~ 100 GeV per nucleon 0, 0, S|- it is known from lattice QCD simulations, for 
, temperatures larger than Tc ^ 170 MeV, confined nuclear matter shall build a phase of deconfined plasma of 
quarks and gluons - the QGP phase Q]- Although the collision energies at RHIC are sufficient to build this 
phase, the main goals of the physics of heavy ion collision (HIC) are far more than only to create the QGP. 
. ^ [ One would like to discover the deconfined QGP under thermal and chemical equilibrium in order to study its 
' static properties such as equation of state, temperature, transport coefficients, etc. But, since it is known that 
^ I the system at RHIC evolves within time duration of the order of 10 — 100 fm/c, it is necessary to study at 
- - ' the same time its dynamical properties. Filling the large gap between the static and the dynamical aspects of 
HIC, relativistic hydrodynamics is one of the most effective methods able to describe the time evolution of the 
expanding QGP Q, and to explore the behavior of nuclear matter in the vicinity of phase transition point 
In order for hydrodynamics to be applicable, the duration of a heavy ion event, r, has to be large compared to 
the equilibration time. It is widely believed that the QGP produced at RHIC behaves as a nearly perfect fluid: 
The transverse radius of an Au nucleus is approximately 6 fm and on the order of 7000 particles are produced 
overall 0]. The motion of particles is relativistic and the nuclei are Lorentz contracted by a factor of 7 ^ 100. 
The duration of a heavy ion event r ~ 6 fm. As it is shown in Q], there is indeed a fascinating agreement 
between predictions from ideal relativistic hydrodynamic models with the experimental data. Whereas plasma 
instabilities, such as electromagnetic Weibel instabilities in relativistic shocks can be made responsible for the 
possible short time thcrmalization of the QGP at RHIC 0], methods based on gauge-string duality starting 
from collisions of gravitational shock waves [§] can describe the enormous entropy production in the heavy ion 
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experiments 

Apart from the confinement-deconfinement phase transition, the chiral symmetry restoration plays an impor- 
tant role in the QCD phase transition. Recent lattice QCD results [10] predict a critical temperature ~ 180 — 200 
MeV for the transition from a chirally broken to a chirally symmetric phase, the chiral phase transition. At 
this critical temperature, following a first order phase transition, the energy density as a function of temper- 
ature is shown to increase suddenly by ~ 1 GeV/fm^. The same phenomenon is believed to occur during the 



QCD phase transition in the Early Universe Here, the phase transition, being of first order, is usually 



followed by a mechanism of supercooling 11| . Depending on the degree of supercooling the kinetics of domain 
formation and growth can be described by the mechanism of nucleation or spinodal decomposition. In contrast 
to the Early Universe, where the QCD phase transition was driven by nucleation of bubbles of true vacuum 
inside the metastable phase, the QGP created in RHIC expands several orders of magnitude faster than the 
primordial Universe, and, bypassing the nucleation process enters the domain of spinodal decomposition ^1^. 
The phenomenon of sudden hadronization is also suggested by several results from CERN-SPS and BNL-RHIC 

The mechanism responsible for the onset of instabilities in an expanding QGP at nonzero temperature 
and baryon chemical potential is studied in [l^ . Using the effective potential of the linear cr-model as the 
thermodynamic (effective) potential, a chiral hydrodynamic description of the expanding QGP near the chiral 
phase transition is derived (see also [l3|)- The effective potential is determined by integrating out the quark 
degrees of freedom that play the role of a heat bath for the chiral fields a and tt with temperature T and baryon 
chemical potential fj.. Instead of using the thermodynamic potential to calculate the total pressure and energy 
density and to obtain the conserved energy-momentum tensor for the expanding perfect fluid, the authors in [12j 
adopt the variational formulation to obtain the hydrodynamic equations for their system 15]. This approach 
provides a natural way to merge chiral and fluid dynamics. Performing a stability analysis up to flrst order, the 
dispersion relation for a plane wave propagating in the above QGP coupled to chiral fields is derived. Solving 
the dispersion relation, the pressure (sound) modes and chiral modes are further determined. 

In the present paper, we will use the same variational method to study the effect of large magnetic fields on 
the sound propagation of a plane wave propagating in an expanding magnetized QGP coupled to chiral fields 
at finite (T, /i) and in a strong magnetic field. For a magnetic field aligned in the third direction, we will, in 
particular, determine the anisotropy in the velocity of a plane wave propagating at angles {(p, 9) of spherical 
coordinates from the external magnetic field. Similar frequency dependent anisotropy was previously observed 
in magnetic fluids in condensed matter physics ^^.^ To mimic the hot and dense QGP in a strongmagnetic 
field, we will use the one-loop effective action of the Nambu-Jona Lasinio (NJL) model of QCD [17|, ll8| at 
finite (T, A*) and in a strong magnetic field B. As it is shown in [19], the one- loop effective action of the NJL 
model exhibits a dynamical chiral symmetry breaking due to the well-known phenomenon of magnetic catalysis. 
According to this phenomenon, in the limit of strong magnetic fields, NJL dynamics is dominated by the lowest 
Landau level (LLL), where even at the weakest attractive interaction between fermions the chiral symmetry of 
the theory is broken by a dynamically generated fermion mass. This mass is shown to depend on (T, /i) and 
the strength of the applied magnetic field. 

The magnetic catalysis has applications in both cosmology [20| and condensed matter physics [2H , and 
seems to be also relevant in the physics of heavy ion collision: As it is reported in Bin, Q, in off-central 
collisions, heavy ions possess a very large relative angular momentum and create very strong magnetic fields. In 
this situation, the presence of topological charge was predicted to induce the charge separation with respect to 
the reaction plane. Large magnetic flelds play therefore an important role in the physics of non-central heavy 
ion collisions and provide a possible signature of the presence of CP-odd domains in the presumably formed 
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QGP phase 13, 23, [24| . The detailed theory of this chiral magnetic effect describing the interplay between the 
chiral charge and the background magnetic field has been developed in [2^.^ Relativistic shock waves, that are 



^ Magnetic fluid is a colloid of tiny (lOOA) magnetic particles or grains suspended in a carrier fluid such as water. The magnetization 
of the fluid varies with the applied magnetic field, typically reaching a saturation of 10'^ — 10^ Gaufi [T^ . 

In [25l an experimental observable sensitive to the chiral magnetic effect has been proposed and the flrst preliminary results have 
been reported by STAR Collaborations 12^ . As it is suggested in [23l . the generation of chirality in the QGP, responsible for 
the charge separation, is the QCD counterpart of the generation of baryon asymmetry in the electroweak phase transition in the 
Early Universe [27 l[2^ . 
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believed to be build in the heavy ion collisions, can be viewed as a possible origin of the magnetic field generation 
in the heavy ion collision. The mechanism is well-known from astrophysics, where the generation of large 
magnetic fields in relativistic shocks plays an important role in the fire-ball model for Gamma-ray Bursts [29[. 
This model proposes that the non-thermal radiation observed in the prompt and afterglow emission from the 
Gamma-ray Bursts is synchrotron radiation from coUisionless relativistic shocks. In collisionless shocks, plasma 
instabilities can generate magnetic fields. Within the shock transition layer the relative motion of the mixing 
pre- and post-shock plasma produces very anisotropic velocity distributions for all particle species concerned. 
Fluctuating electromagnetic fields deflect the incoming charged particles and act as the effective coUisional 
process needed to complete the shock transition. These fluctuating fields occur naturally because anisotropic 
velocity distributions are unstable against several plasma instabilities, such as electromagnetic Weibel instability 
(3o| . The latter is an instability of the currents that result from charge bunching in the beams, and leads to 
spontaneously growing transverse waves with |B| > |E|. In a relativistic shock, where the velocity of the pre- 
and post-shock plasma approaches the velocity of light, the Weibel instability dominates, because it has the 
largest growth rate. In heavy ion collisions, Weibel instabilities are made responsible for the short equilibration 
time of the system from the initial colliding stage 3l|.^ 

The effect of a strong magnetic field to modify the nature of the chiral phase transition in QCD is recently 
studied in H, 34 , S^] . In 33 1 , using the one-loop effective potential of the linear cr-model, the hot and dense 



QGP is simulated and it is shown that for high enough magnetic fields, comparable to the ones expected to 
be created in non-central high energy heavy ion collisions at RHIC, the original crossover is turned into a first 
order transition. In 3J] , the finite-temperature effective potential of linear cr-model in the presence of constant 
and weak magnetic field, including the contribution of the pion ring diagrams in the background of constant 
classical cr-fields is calculated. It is shown that there is a region of the parameter space where the effect of ring 
diagrams is to preclude the phase transition from happening. Inclusion of magnetic field has small effects and 
becomes more important as the system evolves to the lowest temperatures.^ In [s^, the response of the QCD 
vacuum to an external Abelian chromomagnetic field in the framework of a nonlocal NJL model with Polyakov 
loop is studied and a linear relationship between the deconfinement temperature and the squared root of the 
applied magnetic field is found. 

The next question is how the magnetic field would affect the hydrodynamical quantities of an expanding 
magnetized QGP. One of these quantities is the velocity of sound modes of a propagating plane wave in this 
medium. We investigate this question in the present paper. Further, we are interested on the hydrodynamical 
signatures of a chiral phase transition once the temperature reaches the critical temperature Tc- Previous studies 
on the temperature dependence of sound velocity Vg shows that at T = the sound waves has minimum velocity 



37| . As it will be shown in Sect. IVIl at finite (T, /i), apart from an anisotropy in the sound velocity with respect 



to the direction of the magnetic field, the same behavior occurs for a plane wave propagating in the transverse 
direction to the external magnetic field. In contrast, for the waves propagating in the longitudinal direction to 
the magnetic field, the sound velocity has a minimum for T < T^, reaches its maximum at T ^ and remains 
constant after the temperature passes the chiral critical point, i.e. for T > Tc- On the other hand, at nonzero 
baryon density, the sound waves seems to die out, as Vg has a real and an imaginary part. Whereas fHe(ws) 
has the same behavior as Vg for ^ = 0, the 3m{vs) is several orders of magnitude smaller than d\t{vs) and 
oscillates. This behavior which can be interpreted as the onset of the aforementioned spinodal effects is also 
observed recently in j38j , where the dynamical density fluctuations are studied around the QCD critical point 
using dissipative relativistic fluid dynamics with no chiral fields included. This is interpreted experimentally as 
a possible fate of Mach cone at the chiral critical point. The disappearance or suppression of the Mach cone 
would be a signal that the created matter has passed through the critical region, showing the existence of the 
QCD critical point [sst . 

The paper is organized as follows: In the first part of the paper (Sects. HIl and HIT)) . we will derive the 
effective Lagrangian density of the NJL model, consisting of the effective kinetic and the one-loop effective 
potential of the theory at finite {T,fi) and in the presence of strong magnetic field B. This part can be viewed 



^ Recently, an alternative scenario based on Nielsen-Olesen instability 32] which is characteristic for the configuration of a uniform 

magnetic field is introduced in Q. 
* The effect of ring diagrams on the dynamical chiral phase transition of QED in the presence of strong magnetic field is studied 

in H. 
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as a generalization of the results presented in 1^ 39 1 and the methods introduced in [i^ and 41 1 to the case of 
nonzero chemical potential. The effective potential and effective kinetic terms are determined in Sects. Ill Al and 
III CI respectively. The effective kinetic term is in particular determined using a derivative expansion following 
the standard effective field theory methods presented in (4^ . Solving the corresponding gap equation that arises 
from one- loop effective potential of Sect. Ill Al the dynamical mass of the NJL model generated in the presence 
of strong magnetic fields at finite (T, /x) is calculated in Sect. Ill Bl Note that the dynamical mass, being the 
configuration that minimizes the effective potential, can be viewed as the equilibrium configuration, once the 
instabilities in the magnetized QGP are set on. Comparing to the effective action of the linear tr-model used in 



I4I , both the effective kinetic term and the minimum of the effective potential of the NJL model in a strong 
magnetic field depend on (T, fj,) and the constant magnetic field. To have a link to hydrodynamics, we will 
determine in Sect. Illll the energy- momentum tensor T^"^ of the effective NJL model in the presence of a strong 
magnetic field using the effective Lagrangian density from Sect. |lll Here, we will introduce a p olarization 



tensor ilf^" containing the magnetization M of the medium. The same tensor appears also in [43|, |4j| , where 
a MHD description of an expanding dissipative fluid is introduced to study the Nernst effect in the vicinity of 
superfluid-insulator transition of condensed matter physics.^ 

In the second part of the paper (Sects. IIVIIVI|) . we will use the method introduced in [1^] and the results 
from the first part of the paper to present a chiral magnetohydrodynamic description of a magnetized perfect 
(non-dissipative) QGP coupled to chiral fields of the effective NJL model. This is in contrast to the work done 



in 4J] and [4a | where in the presence of a weak magnetic field, the magnetized fluid does not involve any chiral 
fields. In Sect. IIV Al we will first generalize the thermodynamic relations to the case of nonzero magnetic 
field. In Sect. IIVBI we describe the variational method used in 12] by comparing first the polarization tensor 
M'^'^ of the NJL model with the expected magnetization of the magnetized medium (Sect. IIV B ip . and then 
the energy-momentum tensor of the effective NJL model with the corresponding T^'^ of an expanding QGP 
coupled to chiral fields (Sect. IIVB2p . In Sect. |Vl using the hydrodynamical equations, the energy-momentum 
conservation relation and the conservation relations for the number and entropy densities, we will perform a 
first order stability analysis and derive the corresponding dispersion relation of the magnetized QGP coupled 
to chiral fields a and tt. The sound velocity is then derived in Sect. IV Bl We will use the method in [l6| to 
determine the anisotropy in the sound velocity. In our specific model studied in this paper, a wave propagating 
in a plane transverse to the external magnetic field, the anisotropy is independent of the angle if of the spherical 
coordinate system. On the other hand, a ^-dependence arises in the anisotropy function, when the propagating 
wave is in the longitudinal plane with respect to the direction of the magnetic field. In Sect. IVII performing 
a numerical analysis, we will visualize our results in several figures. Our results are summarized in Sect. IVIII 
In Appendix [XI we will generalize the Bessel- function identities appearing in the high temperature expansion 



of Matsubara sums of finite temperature field theory presented in [40| to the case of finite chemical potential 



In Appendix [B] a generalization of the Mellin transformation [4l| for the case of nonzero chemical potential is 



presented. A detailed derivation of thermodynamic relations used in Sect. |V]is presented in Appendix [Cl 

II. NJL MODEL AT FINITE T, /j, AND IN A STRONG MAGNETIC FIELD 

We start with the action of the NJL model in 3 + 1 dimensions at zero temperature and zero baryonic density 

£ = i [v;, {irD^)i^] + 1 [ (H^y + i^pii'i^)'] ~ T, (II. 1) 

where = \{F^vY is the gauge kinetic term. Defining the covariant derivative by = — ieA'^^, with 
the gauge field in the symmetric gauge A'^* = | {0,—Bx2,Bxi,0), a constant magnetic field aligned in the 
xa-direction is generated. In the above symmetric gauge, the only non-vanishing elements of the field-strength 
tensor are F12 = -F21 = B and F^^, = F'"". We have therefore = Note, that (jlLip is chirally 



invariant under C/l(1) x Ur,{1) group transformation. As it is described in 19|, the above theory is equivalent 



^ In condensed matter physics, the Nernst coefficient measures the transverse voltage arising in response to an appUed thermal 
gradient in the presence of a magnetic field |44| . 



4 



to the bosonized Lagrangian density 

£ = 1 [Vi, {il^D^)ii\ + i-f^n) ^-^{a^+ n^) - (II.2) 

where the Euler-Lagrange equations for the auxihary fields a and tt take the form of constraint equations 

cr = -G {^^) , and TT = -G {^i-f^ij:) . (II.3) 

Integrating out the fermion degrees of freedom, we obtain the effective action for the composite fields p = (cr, tt) 
at finite temperature T and density and in the presence of a constant magnetic field eB.^ It is given by 

r[p; eB, ^i,T]^ J d^x C^„{p; eB, T,fi)^ J d^x {Ck ~ fl ^ J^) . (II.4) 

In the following, we will first introduce the effective potential Q, and then determine the kinetic term Ck using 
the method introduced in [l^. 

A. The effective potential of the NJL model at finite T, n in a strong magnetic field 

The effective potential of the NJL model consists of a tree level part and a one- loop part. The latter is determined 



in 19| at zero temperature and density, and in the presence of a constant magnetic field. It is expressed through 
the path integral over fermions 

exp (it{p; eB)^ = J DtpVip exp Q [i), {ij^D^ - cr - i757r)'(/'] ^ 

= exp (Trln {i-f'^Df, -a- . (II.5) 

Using the definition of Tr(- • • ), F can be given by 

r ^ J d'^x V^^'>{p;eB,T = p = 0), where V^^^ = -itT\n{ijf'D^ - cr - ij^ir) . {11.6) 

Here, tr(- • • ) is to be built only over the internal degrees of freedom. It can be further evaluated using the 
Schwinger's proper-time formalism 19|, |46|. The effective potential of the NJL model at zero T and p, fl 



including the tree level contribution V^^\ and one-loop effective potential V^^\ is therefore given by 

n{p; eB,T = p = 0) = V'-°^ (p; eB,T = p = 0) + (p; eB,T = p = 0), (II.7) 



where 



1/(0) = -g-, and = |4 r %coih[eBs)e-^P\ (II., 



Here, p'^ = cr^ + tt^. In the limit of very large magnetic field eB oo,"^ the UV cutoff A can be replaced by 
A ^ As = V eB, and coth(ei3s) ~ 1. The integral over s in (|II.8|) can therefore be performed using 

/>oo 

T{n,z) = dt t"-^e-K (H.9) 



^ It can be shown that the presence of a constant magnetic field, defines a natural scale £g = } , where e is the electromagnetic 
coupling constant. 

^ The strong magnetic field limit is characterized by a comparison between eB and the momenta of the particles included in the 
theory. In this limit |k|||,|kx| <^ VeB, where, |k|| |anii|kj^ | are the longitudinal and transverse momenta with respect to the 
direction of the constant magnetic field, respectively. 
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It leads to 



2G ' 87r2 8^2 (1 In 



O 



eB ' 



(11.10) 



where the asymptotic expansion of r(— 1, z) fa i — 1 + 7^; + In z + 0(z) for z — > 00 is used. Here, 7^; ~ 0.577 is 
the Euler-Mascheroni number. 

Similar methods can be used to determine the effective potential of the NJL model at finite temperature 
and density, and in the presence of a constant magnetic field 47, 4^ 



00 r\ ( 27r fis I 



coth [seB) e 



Here, (3 is the inverse temperature, /3 = ^, and 

oc 

e2(w|r) = 2 ^ e'--("+^)' cos ((2n + 1) u) , 



(11.11) 



(11.12) 



is the elliptic 0-function of second kind. Using the identity [47 1 



where 



e^ulr) = 1 + 2 ^(-1)" e^™'^ cos {2nu) , 



(11.13) 



(11.14) 



is the fourth Jacobian 0-function, and setting coth(ei3s) « 1 in the strong magnetic field limit eB -^00,^ the 
one-loop effective potential in (jH.lip can be separated into two parts. 



{p; eB, T, /i) = v'^^^ (p; eB.T = p = 0) + V'^^^ (p; eS, T 7^ 0, /i 7^ 0). 



(11.15) 



The effective potential fl{p;eB,T,p) including the tree level yf*'^ and the one-loop effective potential T^^^^ is 
therefore given by 



n{p;eB,T,p) = 



p2 eB f°° ds _ 



1 + 2 ^(-l)"e-^ cosh(n/x/3) 



n=l 



To evaluate the integration over s, the integral pi.Op and 

oc 



J dx x"-^ exp \^~^ - 7a;j = 2 ' i^^/h) 



(11.16) 



(11.17) 



can be used. The effective potential (IH.16|) is therefore given by 

p^eB 



nip;eB,T,fj.) 



2G ' 8^2 -r( I'^i 



p \ , P^^ Y^^ZlLK^^pep) coth iin(3) . 



(11.18) 



* At finite temperature T, the strong magnetic field Hmit is characterized by <g eB. 
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where p = \p\. The (T, /i) independent part of pi.lSp is exactly the same as pi.lOp and can be similarly 
evaluated in the strong magnetic field limit. As for the (T, ^) dependent part, it can be expanded in the orders 
of (p/3), using the Bessel function identities from Appendix 1X1 In particular, we will use the identity 



E 



i-iY 1 

^--^Xi(fe)cosh(fe') = -z 
I 8 



l-27B-21n - 



(-1)"(1 + n) / z 



§ (^(2 + r^))2 



(2n) 



IZ 

2^ 



2tt) 



fe=i 



22fe+l 



/ / \ 2k -J 

i^j C(2fc + l)+2;{Li2(e^("+*"'^) +Li2(e-*("+'^'))}, (11.19) 



with z = p(3 and z' — jjLjS from (|A.23[) to evaluate the sum over i in (jll.lSp . In pi.l9p . the polygamma 
function il^{z) = ^ Inr(z), and (z) = ^^'^ is the (m + 1)*'' logarithmic derivative of the r(z)-function. 
Furthermore, the dilogarithm function Li2(z) is defined by 



k=l 



and satisfies 



(11.20) 



Li2(z) 



ln(l - t) 



dt. 



(11.21) 



The effective potential of the NJL model pi.l8|) at finite (T, /i) and in the limit of eB ^ oo is therefore given 

by 

ni- urr ieB f p^eB ( f eB(3^ 



eB 



E 



2G 87r2 87r2 
(-1)" (n + l)(p/3)2("+i) 



87r2/32 (167r2)" (r(2 + n)Y 
eB{pPf (-l)'=(22fe+i - 1) (pP 



37r2/32 



E 

fe=i 



22fc+l 



C(2fc + l) + ^{Li2(e'(-- 



(11.22) 



It exhibits a spontaneous symmetry breaking, which is mainly due to a dynamical mass generation in the regime 
of lowest Landau level (LLL) dominance. The type of the phase transition as well as the critical temperature 
depends on the strength of the magnetic field eB. Fig. 1 shows two examples of 57 for a) eB — 10 GeV^ which 
is equivalent to B '--^ lO'^^ Gaufi,^ and for b) eB = 10""* GeV^ which is equivalent to B 10^^ GauB. In both 
cases the potential 51 from pi.22[) is calculated up to 0{{pj3y) and 0{{p(3y), and p(3 is chosen tohe p — 10~^ 
GeV. Let us note that the qualitative picture that arises here does not depend too much on p. 

In the next section, the effective potential pi.22|) will be used to determine the dynamical mass, that is 
generated due to the phenomenon of magnetic catalysis in the strong magnetic field limit 19(|. 



B. Dynamical mass of the NJL model at finite T, /i in a strong magnetic field 

As it is shown in [l9| . a constant magnetic field in 3 + 1 dimensions is a strong catalyst of dynamical chiral 
symmetry breaking. This phenomenon leads to the generation of a fermion dynamical mass even at the weakest 



As it is shown in (361, eB GeV^ is equivalent to B = 1.691 X 10^" GauB. 
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a) n(p/T) for cB = 100 GeV^ and /i/T = 10"' 



a) n(p/T) for eB = lO"'' GeV^ and ^i/T = 10"^ 




FIG. 1: a) Effective potential Q. of the NJL model as a function of p/T for eB = 10 GeV^ and = 10~^. Here, a 
second order chiral phase transition occurs at Tc ~ 100 MeV (far right), b) Effective potential Q of the NJL model as 
a function of p/T for eB = 10~* GeV^ and fiP = 10~^. Here, a smooth crossover occurs at Tc — 100 MeV (far right). 
Here, the NJL coupling G = 5.86 GeV~^ and the sum over n in (|H.22p is up to n, fe = 7. 



attractive interaction between the fermions. The essence of this effect is a dimensional reduction from D — > D — 2 
in the dynamics of fermion pairing in a magnetic field. At zero temperature and vanishing chemical potential, 
the dynamical mass can be determined using the effective potential 



from (jlLlOp through the gap equation 

dn{p;eB,T = ^ = 0) 



da 



-4 



0. 



(IL23) 



(To#0,7ro=0 



Here, the minimum of the potential is assumed to be at unknown but constant cto and vanishing ttq. The gap 
equation pi.23|) leads to 



G ~ 47r2 



r 0, 



eB 



(11.24) 



that, after excluding the trivial solution (Tq = 0, and using the approximation lim r(0, z) ~ —7b — Inz, leads 
to the non-vanishing dynamical mass at zero temperature and chemical potential 

^'^ \ , with Cm ^ e"''^/2 ~ 0.74. 



ao{eB, T = /i = 0) = Cm v eB exp 



GeB 



(11.25) 



The same result arises also in |39j . Next, we will determine the mass gap at finite (T, /i) and in the strong 
magnetic field limit. Using f2 (p; eB, T, pi) from (jn.l6[) . we arrive first at the gap equation 



= 



dn{p;eB,T,p) 



da 



o"o5^0,7ro=0 



G 47r2 



e=i 



(11.26) 



The temperature independent part of pi .26p can be evaluated as in pi.24p . leading to the iterative solution 

/ A •} ' ' ^ 



ao{eB,T,fi) = CmV eB exp 



GeB 



(n.27) 



where the numerical factor Cm — 0.74 is in pi.25|) . Setting p, = 0, the above result pL27|) coincides with the 
dynamical mass determined in [i^. At this stage, we would however use the Bessel- function identity (|A.16p 
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from Appendix [XI that includes a complete expansion in the orders of {(Jof3) and 

2k 



fc=l 



, V (-1)" ^^0/3 

n=l ^ 



(2n) ( 1 

2 2 



TT 



2") (l + ^ai. 

2 2 



(11.28) 



We arrive at the dynamical mass of the NJL model at finite temperature and density and in the strong magnetic 
field limit 



ao{eB,T,n) = C,„VeB exp 



'Tol3\ ^ (-l)^(2^^-+^-l) 



2k 



2^ Zl2n 



•1=1 ^ ' ^ 



(2n) 



(11.29) 



In pi.28p . the polygamma function ij^iz) is defined as in (jll.f 9^ . The identity pi.28[) is a generalization of 
Bessel-function identities from 40 1 for /i 7^ and will be proved in Appendix \K\ 



o-o (T, fi) in GeV for a fixed eB = 1 GeV 
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FIG. 2: Qualitative dependence of NJL dynamical mass cro from (|IL29|| on T and in a strong magnetic field eB = 0.1 
and eB = 1 GeV^. The dynamical mass abruptly decreases by a change of the magnetic field of one order of magnitude. 

Fig. 2 shows qualitatively the (T, /i) dependence of ao{eB,T, fi) from (|II.29|) for eB = 0.1, 1 GeV^, corre- 
sponding to i? ~ 10^^, 10^° Gauf5, respectively. Since the solution (jll.29p is an iterative one, it is necessary to 
choose an initial value for cto on the r.h.s. of pi.29p . In Fig. 2, this initial value is chosen to be ao{eB, T = /i = 0) 
from pi.25p . Further, the sum on the r.h.s. of pi.29p is up to n = fc = 7. The NJL coupling G is chosen to 
be G = 5.86 GeV~^. Fig. 2 shows that the dynamical mass vanishes for a fixed value of eB by increasing 
the temperature. Its qualitative behavior does not depend too much on the chemical potential fi. In the next 
section, the kinetic term of the effective action pi.4p will be determined in a high temperature expansion. 



C. The kinetic term of the NJL model at finite T, /i and in a strong magnetic field 



The kinetic term Ck of the effective action pi.4p at zero temperature and density, and in the presence of 
a constant magnetic field is previously determined in [l9| . In the first part of this section, we will briefly 
outline the arguments in 19| . We then generalize them for the case of non- vanishing temperature and chemical 
potential. 

Let us start with the general structure of the kinetic term of the effective action in a derivative expansion 



(11.30) 



This step turns out to be useful in Sect. |V] 
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that is implied by the Ul{1) x Ufi{l) symmetry of the original Lagrangian density. Here, p — (ct, tt), and 
F^^ ,i = 1,2 are structure constants depending nontrivially on the external magnetic field B. Following the 
general field theory arguments, and the definitions 



5a{x)5a{Q) 



<To#0,-n-o=O 



S'k{x)St:{0) 



<To740,7ro=0 



(11.31) 



the structure functions can be given as 



F^ 



- I d^xxf^x" 



1 



Fr = -4 / d^xx^x" 



S7r{x)STT{0) 
Sa{x)da{0) 2 



1 p^i. 



(11.32) 



where T = — iTr \n{ij^Df^ — cr — ijzir) is the effective action from (jll.Sp . In (jll.3ip the minimum of the potential 
is chosen to be at cto = const, and ttq = 0. To calculate F/"', i — 1,2 from (jll.32[) . we use the definition of the 
fermion propagator iS^^ — i^^Df^ — ctq/^ and arrive therefore at 



5a{x)5(j{Q) 

52f 



(57r(a;)57r(0) 



-itr {S{x,Q)S{<d,x)) 
+itr(^(x, 0)75^(0, x)75). 



(11.33) 



Plugging pi.33p back in (|II.32p . the structure functions in the momentum space are given by 



Fr = +77 



^2 



d*k 

(2^ 

d^k 
(2^ 



tr ^(fc)75 



tr S{k) 



dkp^dky 



75 



dkndki, 



-:F^ 



(11.34) 



As it is shown in 19|, the fermion propagator in a constant magnetic field and at zero {T,fi) is given by 



S{x, y) = exp ( - (x - yf A^^^x + y) ) 5(x - y), 



with the Fourier transform of 5* 



(11.35) 



(11.36) 



where n labels the Landau levels that arise in the presence of a constant magnetic field, and ky as well as kj_ 
are the longitudinal and transverse momenta with respect to the external magnetic field, i.e. for B = Be^ we 
get k|| = (fcg, k^), kj^ = (fci, fc2). Further, Dn(eB, k) is given in terms of Laguerre polynomials L„ by 

D^ieB, k) = (k|| • 711 + ao) |o_i„ [2^^ - 0+L,,^i (^2^) | - 2k^ • 7||4^2i (^2^^ , (11.37) 

where 0± = ^ (1 ± Z7i72sign(ei?)). In flS], the fermion propagator in a constant magnetic field (|II.35|) - pi.36|) 
is used to determine = 1,2 from pi.34p . But, since we are interested only on the strong magnetic field 



See e.g. in [42l | for more details. 

^■^ Note that here, o-q = const, plays the role of the dynamical mass. It will be determined analytically by solving the corresponding 
gap equation in Sect. Ill Bl 
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limit eB — )■ oo, where the fermion dynamics is dominated by the lowest Landau level (LLL) pole of the fcrmion 
propagator, it is enough to use the fermion propagator in the LLL approximation 

^(x)=p^(x^),s||(x||)0_, (11.38) 

with 

PiXi=— e and 5||(x||) = / — J^^- ^. (11.39) 

2tt J {2iiy (k|| • - (Jo) 

It arises from (|II.35P - pi.37p by setting n = 0. In Sect. IIVI a hydrodynamical description of the NJL model in 
the presence of a strong magnetic field will be presented, where the combination G^'' = F-j"'' + 2^^"^ is shown to 
be relevant. Plugging therefore pi.38P - (|II.39p in (jll.34p . the structure function G^'^ in the LLL approximation 
are given by 



= G- = -^/^tr(l„(k„)0_S„(k„)0_), 



where S||(k||) — (k|| • 7|| — ctq) ^ is the Fourier transform of S||(x||). Plugging S||(k||) in G'"' and using the 
identities tr(l) = 4, tr(7"7'') = 4g'^^ ^ 0_,0_7, = 7„0_, tr(0_7,0_76) = 2gab, and tr(0_7,) = for 
a, 6 G {0, 3}, as well as tr(C'_) = 2, we arrive first at 



and 



^eB f <Ph 1 /„,9 „ 9 9. Mm+'^l^ 



''''^--JW)2 (kp^ ^^^^ + (^i + - ^Tkp^ ) ' (^^-^^^ 

whereas the non-diagonal terms identically vanish, i.e. 

TT J (27r)2 (kjj - 0-^)3 \^ (kjj - <7^) J 

To calculate G^'' from pi.40p at finite temperature and density, we will first perform the following necessary 
replacements 

(2£ + l)7r 

fco — > iko with fco — {oje — iji) , and LOg — , (11.44) 



as well as 



/CXj 
(fk\\ -^iT ^ / dfcg, (11.45) 



t=-oo 



and arrive at 



G" = G22 



oo 



dfca A:o + A:| - CTq 



eB ^ f dk3 (fc| - fc^ - cr;! + 6ct§A:2) 



(-,00 ^ 

^33 _ f dks (3fc| - 6fc|ag ~ 4fc|fcg + fcp^ - al) 
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To evaluate the /cs-integration as well as the sum over the Matsubara frequencies we use the Mellin transfor- 
mation technique. This technique, which is originally introduced in [4l'| for finite temperature and zero density, 
is generalized in Appendix |B] for finite temperature and nonzero chemical potential. Using in particular the 
main identity including an expansion in the orders of (jnP), 



1 



E 

— — o 

CX 

k=0 



d^k (k^)°fc^* _ 1 r(|+a) /27r 

(2^)'^ (fc2+k2 + m2)" ~ 2(47r)'^/2r(a)/? T (|) \ (3 



-2a+d+2{t+a) 



"2^ 



2k 



(11.47) 



and neglecting the temperature independent singularities, we arrive at 

\2 r 



1 



327r3 

15 (ao/3) 
40967r7 



1 

4 r 



3 i<Jof3y 



C 5 



1 



"27 



2567r5 



^ ' ^' 2 2^ 



(11.48) 



and 



G' 



00 



2567r^ 
295(ao/j) 
1287r4 



^ ' ^' 2 2^ 



+ (Ai^ - 



'2 27r 



(m 



23(/3(7o)2 



87r2 
f O ((^To/9) 



C 5; 



1^13 



(11.49) 



and 



G33 = 



5127r5 
55(ao/3) 



+ (m->- 



+ 



1287r4 



^ ( 7 1 _ !^ 
'2 2^ 



87r2 

-0((ao/3)6). 



+ (/^ -m) 



(11.50) 



In pi.47P - (III.50|) . C(s; a) is the generalized Riemann zeta function given by C(s; a) — X]fc°=o(^ + where any 
term with (fc + a) = is excluded. In pi.48[) . the digamma function V-'(^) = ^ Inr(z) arises by regularizing 

C(i;a), 



lim C(l + e; a) = lim — 

£^0 ' e-^o e 



V'(a). 



(11.51) 




FIG. 3; Qualitative dependence of = G", i = 0, • ■ ■ , 3 on (T, /i) for fixed eB = 1 GeV^ corresponding to B = 10' 
GauB. They exhibit a qualitative change in the transition region Tc ~ 100 MeV and small /i. 
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Fig. 3 shows G' = G",i = 0, - ■ ■ ,3 m the (T, n) phase space for eS = 1 GeV^ or equivalently B - 10^" Gaufi. 
To determine G% i =, 0, • • • , 3, we have used the expression PL29I) for aQ{eB, T, fi). The iteration is started as 
in the previous section with the value of the dynamical mass at zero temperature and density (jQ^eB, T = fj, = 0) 
from ((TL25)) . The smu in (|IL29|) is taken up to n = fc = 7 and G = 5.86 GeV"^^ To determine G\ i = 0, • • • 3 from 
(|II.48|) - pi.50p . we have used pi.47|) . For the structure functions in Fig. 3, the high temperature approximation 
in pi.47p is up to order (cto/?)^- All three structure functions G',i = 0, • • • , 3 exhibit a qualitative change in 
the vicinity of the critical temperature Te- 
la the next section, we will use the effective kinetic term pi.30|) and the effective potential (III.22[) to determine 
the energy-momentum tensor of an effective N JL model described by the effective action (111.41) . 



III. THE EFFECTIVE NJL MODEL AT FINITE T, fi AND IN A STRONG MAGNETIC FIELD 

The main goal of this paper is to present a hydrodynamical description of an expanding chiral and magnetized 
QGP at finite temperature T and baryon chemical potential fi. It will be described using the NJL model of 
QCD in the presence of a constant magnetic field, that exhibits, as we have shown in the previous section, a 
chiral symmetry breaking due to a dynamically generated fermion mass. The thermodynamics of the NJL model 



has been extensively studied in the past few decades [17|, ll8|. A systematic study concerning the existence, 



location and properties of the critical/tricritical point of various NJL type models is recently presented in [49|. 
In all previous treatments, the effective potential of the theory, il, was of fundamental interest. Within a mean 
field approximation, f2 can indeed be identified with the thermodynamical partition function of the theory, 
Z, through the relation 0((Jo,T, ^) = —-^InZ, where (3 — is the inverse of the temperature and V is 
the volume. Using this partition function, it is then easy to derive all thermodynamic (static) quantities, like 
pressure, entropy, etc (see next section for more details on thermodynamics). 

We, however, are interested in the dynamical properties of an expanding magnetized QGP near the chiral 
critical point. To study the time evolution of such a system, we shall go beyond the mean field approximation 
and to consider the full effective action of the NJL model pi.4|l including the effective potential pL22p as well 



as the kinetic term pL30[) . The latter depends, in contrast to the linear a model used in 1^, on (T, /i) and 



the external magnetic field. The linear a model is in particular used in IJ, Il4| to describe an expanding QGP 
coupled to chiral fields. 

In this section, using the full effective action of the NJL model, we will derive first the equation of motion 
(EoM) and eventually determine the energy-momentum tensor of the effective NJL model at finite (T, /i) and 
in the presence of constant magnetic field B. The energy-momentum tensor will be then used in Sect. IIVI 
to present a magnetohydrodynamical (MHD) description of an expanding perfect QGP near the chiral critical 
point. 

Using the effective action pl.4p . and in particular the general structure of the kinetic term (III.30p . the EoM 
of the chiral fields p= (cr, tt) of the strongly magnetized NJL effective model can be derived from the ordinary 
Euler-Lagrange equation 



It is given by 



where 



apk o{d\pk) 



Fl^d^d^Pk + "^F^ i^p^d^d^p, + d^p,d^p, ~ ^ ip^^^P^) iPjd^Pj) ]=-Rk, (III.2) 



Rk ^ 1^. (111.3) 

opk 

Here, the effective potential fl — 1/'"^ + 1/'^^ includes the tree level potential and the one loop effective 
potential V*^^) presented in pi.22p . In pil.2[) . we have used the property F^'^ = F^'^ g'^'''^ (no summation over 
p, V is considered) and denoted as in the previous section Ff^ = F/", with i = 1, 2 and /i = 0, • • • 3. For F^ = 1 
and F2 = 0, this result is indeed comparable with the results presented in [13,13! where the equation of motion 
of a linear a model at finite (T, p) and with no magnetic field is given by Opk = —Rk- 
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Let us now look at the symmetric energy- momentum tensor corresponding to the effective action pi.4[) , that 
arises from the standard relation 

rir 

where the dynamical fields ^ are given by ^ = {pk, A^^*}. The energy-momentum tensor will then consist of 
two parts: T^^^^^ = T^^^ + Tj^^y Here, the matter part ^'(p), and the gauge part T^'^^ are given by 



7^(p) = (^APfe^.P^c) + ^ef^"^ (p,9aP.) {p.d.Pj) - 9^""^, (in.5) 

with 

^ _2i^f g^^y'"" + F^g^'^g^''', i = 1,2, (III.6) 

and 

T^';^ = -g^-'j^ - 2n|f^ = -g'^"^ - ^ A^^^^, (ni.7) 

with T = K-F^i.)^, respectively. The relation ([TTLT)) can be brought in a more familiar form 

Tr;, - (f^F'-- + g^^'-^F.^F-p) ^ + g^'' (-T , (III.8) 



that appears also in J6J and reduces to the well-known Maxwell tensor 



T^^ = -F^F^" ~ gi'^-Fp^F'^P, (111.9) 



once C^tf — For F^ — 1 and F^ — 0, piI.Sp is comparable with the matter part of the energy-momentum 



tensor T/"^ of the linear cr model 

(p) 



T^;^ = i (-2g''^g^'' + g^''/-) (aAPfc^.Pfe) - 1^(^)5^^ (ULIO) 



The energy-momentum tensor (|IIL10p is used in 12, llJ] to present a hydrodynamical description of the linear 
a model. To derive the energy- momentum conservation relation, we first identify the term f ^qP' in (|III.7p 
with the polarization tensor M^'^, that appears also in [44],^'^ 



fir 

Mt-^ = Ff'^^-P^. (III. 11) 
The total energy-momentum tensor can then be given as 

T^'^^^^ = Tj!" - F\M^'^, (III.12) 

with 

T^" ^Tl"^^- g^'^J^, (III. 13) 

where T|/^^ is from (jlll.sp and —g'^'^T is the first term on the r.h.s. of (|TTL7}. Using the EoM of the effective 

NJL model (jIII.2p and the fact that for constant magnetic field d^J- = 0, it is easy to show that 

a^T,^" - 0, (III.14) 



See next section for more detail on the polarization tensor M^'^. 
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and that the energy-momentum conservation equation reduces to 



'-'fJ--^ {tot) 



(III.15) 



The same relation appears also in 4j|, where a magnetohydrodynamical description of the Nernst effect in 
condensed matter physics is presented. In the next section, we will identify the above energy-momentum tensor 
piI.12[) arising from the effective NJL model in a strong magnetic field with the energy-momentum tensor 
of a perfect chiral and magnetized fluid. In contrast to the description in [44|, the magnetohydrodynamical 
description of this theory will include a chiral symmetry breaking arising from the dynamical mass generation 
in the presence of a strong magnetic field. 



IV. MAGNETOHYDRODYNAMIC DESCRIPTION OF THE EFFECTIVE NJL MODEL AT 

FINITE r, n 

In the first part of this section, we recall standard identities of thermodynamics without matter fields which are 
often used in hydrodynamical models with no background magnetic field. In the second part, generalizing the 



method introduced in 1^ , we will implement the chiral fields in a hydrodynamical description of an expanding 



magnetized QGP including a spontaneous chiral symmetry breaking in the presence of a strong magnetic field. 



A. Thermodynamics in a constant magnetic field; General identities 

Let us start with the identity 

U(V,N,S,M) = -PV + TS + ^N + B-M, (IV.l) 

that defines the internal energy, which is an extensive function of extensive variables; the volume V , the entropy 
iS", the baryon number N, and the magnetization vector M. Here, P is the pressure, /i is the chemical potential 
and B is the constant induced magnetic field. Note that in nonrelativistic systems, N is generally the number 
of particles, which is conserved. In a relativistic system, the number of particles are not conserved: it is always 
possible to create a particle-antiparticle pair, provided the energy is available. In this case N no longer denotes 
a number of particles, but a conserved quantity, such as the baryon number (see [l| for a recent review on 
relativistic hydrodynamics). In what follows, we will assume that the magnetization vector is parallel to the 
external B field, i.e. in our setup B = Be^, we get B • M = BM. The differential of internal energy is given by 
the thermodynamic identity 

dU = -Pdy + TdS + ^idN + BdM. (IV.2) 

The identity pV.ip can be derived by differentiating the relation 

U(AV, AN, AS, AM) = AU(V, N, S, M), (IV.3) 

with respect to A, taking A = 1, and using (irvT2|l . Differentiating (|rVT|) and using again pV.2p . we obtain the 
Gibbs-Duhem relation 

ydP = SdT + fidN + MdB, (IV.4) 

that reduces to 

dP = sdT + nd^i + MdB, (IV. 5) 

where s = ^, n = ^, and M = the entropy density, the baryon number density and the magnetization 
density are introduced. Defining further the energy density e = y, pV.ip and (|IV.2p are given by 

e + P = Ts + fin + BM, (IV. 6) 

and 

de = Tds + fidn + BdM. (IV.7) 
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Using the Gibbs-Duhcm relation pV.5p . the thermodynamical variables are given by 

Similarly the temperature T, the chemical potential /i and the magnetic field B can be defined by PV.7P using 

-(S).. -(^). 

In the next section, we will generalize the above thermodynamic relations to the case where the energy density 
e depends explicitly on the chiral fields p — (a, tt) . 

B. Chiral magnetohydrodynamics in and out of thermal equilibrium 

Extended thermodynamic identities 



Following the variational method presented in [I4I , we treat the gas of quarks for the chiral fields a and tt as a 
heat bath with temperature T and baryon chemical potential fi. In thermal equilibrium, this identification is 
justified by the relation 

Po = -nip;eB,T,^i)\^„^o..o=o, (IV.IO) 
which is often used in the standard thermal field theory. Here, Pq is the thermodynamic pressure given by 

Po = -^lnZ, (IV.ll) 

where /? = is the inverse temperature, V the volume and Z the thermodynamic partition function. The 
effective potential il in (IIV.10[) arises by integrating out the fermions from the original theory. The configuration 
(ao ^ 0, ttq = 0) builds the stationary point of the effective potential or equivalently describes the thermody- 
namic equilibrium. In Sect. |Vl where the sound velocity of a plane wave propagating in the magnetized QGP 
will be calculated, we will assume that the same relation between the pressure and the effective action as in 
PV.10|) . i.e. P ^ — eB, T, p), still holds in a system which is out of equilibrium. Here, p = (a, tt) is the out 
of equilibrium configuration of the chiral fields and depends in particular on the space-time point x. 

As we have seen in the previous section, the thermodynamic pressure is related to the energy density e 
through the relation pV.Gp .^"^ 

P = -e + Ts + pn + BM. 

In the above merging process of thermodynamics and quantum effective field theory, expressed by pV.10[) or 
more generally by P = — eB, T, p), it is therefore natural that the energy density e(p; n, s, M) depends not 
only on (n, s, M), as in the ordinary thermodynamics, but also on the chiral field p = (cr, tt). For e = e(p; n, s, M), 
PV.7P is therefore generalized as [I4I 

de = Tds + pdn + BdM + Rkdpk, (IV.12) 

with 

\9sJn,M,P. \dnJ,^M,p, V5^^/„,s,p, \dpkJ„^,^M 



In the stability analysis, that will be performed in Sect. |Vl we will assume that B, T, ^ and M are always constant. We will vary 
only pk = (cr, 7r) , n and s and consequently P and e around their equilibrium configurations ao,no,so as well as Po and to ■ 
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In thermal equilibrium, where the only relevant configuration is (ctq 7^ 0,7ro = 0), we have in particular, 
eo = e(cro; no, so, In this case, the last relation in (jlV.lSp is modified as 



R,,= {-^] , and R^,= {-^] . (IV.14) 

Using further pV.6[) and pV.12p . the corresponding Gibbs-Duhem relation pV.5[) becomes 

dP ^ sdT + ndfi + MdB - Rkdpk- (IV.15) 

This implies the definitions 

^^K-B.p, K^^^Jt.b^p, \dBjT,^^p^ \dpkJT,p,B 

In thermal equilibrium, Po = P{i^q\ eB, T, fi), the last relation in pV.16P is modified, upon using (lIV.lOp . as 

^^^^(^niao■eB,T,,)^ ^ /gl^^; eB, T, ,) X ^ ^^^^^^^ 



Note that the definitions PV.17P are indeed equivalent with pil.3p for a = ao and tt — ttq = 0. Using the 
EoM pil.2p for this constant field configuration, i?o-o f^^d Rj^^ from PV.ITP as well as (IIV.14P vanish, i.e. 



MHD description of the magnetized QGP 

Using the above generalized thermodynamic relations including the chiral matter fields in thermal equilibrium, 
it is indeed possible to derive all the thermodynamic (static) quantities of the effective NJL model in the 
presence of a strong magnetic field. It is, however, the goal of this paper to go beyond the thermal equilibrium 
to study dynamical properties of a magnetized QGP near the chiral critical point. To this purpose, we will use 
relativistic hydrodynamics which can be used as an effective description of the theory in the same footing as all 
the other effective descriptions including the time evolution of the fluid. There is only one assumption, which is 
associated with hydrodynamics: the local thermal equilibrium (LTE). No other assumption is made concerning 
the nature of the particles and fields included in the fluid. Their interactions and the classical/quantum nature 
of the phenomena involved are fully encoded in the thermodynamic properties, i.e. in the equation of state (see 
[l[ for a recent review on relativistic hydrodynamics). 

The fundamental ingredients of a hydrodynamic analysis at LTE are the conserved quantities and their 
equations of motion. Here, generalizing a similar treatment from [44j , where a similar MHD theory in the 
vicinity of superfluid-insulator transition in two spatial dimensions is evaluated, we introduce the four velocity 

= that represents the velocity of the system in LTE with respect to the laboratory frame. In the rest 
frame — (1,0,0,0). It satisfies the normalization u^u^ = 1. Further, we will use the standard definition of 
the metric tensor g^'^ — diag (1, —1, —1, —1). 



Let us start with the conservation laws of total energy-momentum tensor 44 1 

u.d,T^;'^^^^u,F^''4""\ (1V.18) 

where the hydrodynamic energy-momentum tensor in the presence of a background magnetic field is defined by 

r^';;;^) ee (e + P)u^m'' - P.g^'^ + T^" - F\M''^. (IV.19) 



Here, T^^ is the energy magnetization current. It appears also in 4j|, where the authors do not determine it 
explicitly, and use only its property d^Tj^'^ ~ 0. The latter will be shown to be only valid in the fleld free case 
and for P = 0. In our case, however, where, the perfect magnetized fluid is coupled to chiral flelds a and tt, the 
combination T^"' = (e + P)u^'u'' - Pg^"' + Tj^" satisfies ^f,T^"' = 0. 

On the r.h.s. of pV.lSp . the total bary on current density J^^^^j is defined by 

J(';^,) =nu'^-H9.M^", (IV.20) 
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and satisfies the total current conservation law 



(IV.21) 



In pV.20[) . n is the baryon number density, as is defined in the previous section/^ and M^" is the polarization 
tensor (ii-liil 



/o 
















Ms 








-Ms 








\o 








0/ 



Note that due to the antisymmetry of M'^'^, (jIV.21|l implies the third hydrodynamic relation 

d^inu^') = 0, 

which is then completed with the fourth hydrodynamic relation 



(IV.22) 



(IV.23) 



(IV.24) 



The above entropy current density conservation (IIV.24|) is valid only in a perfect (non-dissipative) fluid. Next, 
we have to find a link between the hydrodynamic quantities defined in the above relations and the field theory 
identities derived in Sect. HIl as well as the thermodynamics relations from Sect. IIVBI In what follows, we will 
first define Mg in terms of magnetization density M that appears in the thermodynamic relations, PV.5[) . pV.12[) 
and PV.15[) . This will be done by identifying the polarization tensor Af^^ with the tensor defined in pil.ll[) 
using the standard field theoretical methods. We will then compare the hydrodynamic energy-momentum tensor 
PV.19P ; with the total energy-momentum tensor of an effective NJL model in a strong magnetic field presented 



m P11.12p . and dete rmine T^" explicitly. 



The polarization tensor 



Let us start by understanding the relation between Ms appearing in pV.22p with the magnetization density M 
appearing in the thermodynamic relations PV.5[) . pV.12[) and pv.isp . To this purpose, we define a Gibbs free 
energy density Q in the presence of a constant magnetic field B 5Qi] 



g{p;eB,T,^i) = — 



n{p;eB,T,^,)-HB, 



(IV.25) 



where fl is the free energy density, and H is the external magnetic field. Whereas in vacuum H = B, m a medium 
with finite magnetization density, the external magnetic field H is different from the induced magnetic field B 
50|. Using thermodynamical argument, we will prove the relation B = M + H, that, apart from a normalization 
factor, appears also in [s^]- To do this, let us evaluate Q at its stationary point with respect to p= {a, tt) and B. 
At this point, Q describes all the thermodynamic properties of the system at thermal equilibrium. The stationary 
point, which was in the previous sections characterized by the configuration (cr, tt) ((Tq, 0), corresponds to the 
solution of the gap equations 



dg 

da 

m 

dB 



= 0, 



= 0, 



or 



or 



'dB 



o"o#0,7ro=0 



0, 



H + B = 0. 



(IV.26) 



The first equation in pV.26p is the same as pi.26p . In the second equation, we set, Pq 

dPo 
dB 

B = H + M, 



PV.IOP . Using now -|§|^„5^o,7ro=o = ^ = M from pV.16p . we arrive at 



from 



(IV.27) 



This is in contrast to |44j . where J'f^^^\ is the electrical current and n is the charge density. 
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as was claimed. Comparing now the expression on the r.h.s. of pV.25p with the definition of the effective action 
from pi.4p . we get 



g{p; eB, T, /i) = + Ck - HB. 

At the stationary point, the effective kinetic term vanishes, i.e. Ck\a-o^o,7ro=o 
|(To5^o,7ro=o — from pV.26|) leads therefore to 



dB 



-H = M -B. 



(IV.28) 

0. The gap equation 
(1V.29) 



Using the above relations, we will find in what follows the relation between Alg from (IIV.22|) and the magneti- 
zation density M from pV.29[) . To do this, we will identify the polarization tensor M'^" from pV.22p with the 
polarization tensor defined in the field theory treatment of the effective NJL model (jllLlip , i.e. M^^ = F^^^ ■ 
For B = Be^, we arrive at 










M, 

-Ms 

/ 



0\ /O 0\ 
B 



-B 

V / 



(1V.30) 



that upon using T = \B^ ^ leads to ~ "^^^ . At the stationary point, this relation can be compared with 
PV.29[) and leads to the definition of Ms in terms of the magnetization M and the induced magnetic field B 



Afs|CTo#0,7ro=0 — M — B. 



(1V.31) 



Note that without considering the gauge kinetic term —T in the effective action pl.4p . Ms|cro5^o,7ro=o would be 
the magnetization M at thermal equilibrium, as is also introduced in [44]. 



2. The total energy-momentum tensor T'^'^^jj 

Our next task will be to compare the energy-momentum tensor of the effective NJL model in the presence of 
a strong magnetic field piI.12[) with the energy-momentum tensor of a magnetized fluid consisting of chiral 
fields pV.19p . First, let us show that the energy- momentum conservation law, pil.isp . on the field theory side 
reduces to (IIV.18I) in the hydrodynamics side. To do this, we consider (IIII.12p . that satisfies 

"-5M^^(^ot) = -u.F\M'^\ 

= u,F\J^,^,^, (1V.32) 

as expected. To derive pV.32p . we have used the facts from previous section, and replaced M^'^ from pil.lip ap- 
pearing in the energy-momentum tensor in the field theory side, with M^'^ that appears in the definition of J^^^jj 
from pV.20p on the hydrodynamics side. Further, df^Tj^'^ = with T^"^ from pil.l3p . and the antisymmetry of 
the field strength tensor F^'^ are used. 

Next, we will determine T^", that appears in pV.19p . To do this, let us use the thermodynamic rela- 
tions pV.ep . pV.12p and pv.isp . and the hydrodynamic equations pV.23p and pV.24p to determine first 
Uiydfj, {wu^u^), where w = e -\- P , 

= df, {fin + Ts + BM) + (Ts + ^in + BM) d^.u'' 

= u^'d^P + u^'Rkdf.pk+u.df.iBMu^'u''). (1V.33) 

On the first line, we have used Ui,u'^ — 1, and set Uudf^u'^ = 0. The expression on the last line arises by making 
use of pV.23p . pV.24p . and pV.lSp . as well as the fact that for constant magnetic field d^B = 0. The relation 

u.d^ {{e + P- BM)u^u'') = u^d^,P + uf^Rkd^pk, (1V.34) 
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with zero magnetic field appears also in [l^. It corresponds to the Euler's equation of energy-momentum 
conservation in the non-relativistic hydrodynamics, and will be used in Sect. |Vl to derive the dispersion 
relation and the sound velocity in the vicinity of chiral critical point. Plugging this relation back in (jlV.lQp . 
and using the energy- momentum conservation pV.18[) . we arrive first at 



M.a^rr = -u^Rkd^Pk - u,d^ {BMu^^u") . (IV.35) 

9pk 



Using now the definition of Rk = 4^ from pil.3[) and assuming constant T, fi and B, we then get 



T^" ^ Pg^'^' - BMu'^u'' + a'", (IV.36) 

where P = —ft is used. Here, the constant tensor C^" satisfies Uiydf^C'^'^ = 0. Plugging (jIV.36[) back in (jIV.19[) . 
the total energy-momentum tensor is given by 

'^(tot) ^i^ + P- BM)uf'u'' - Fj^M^^ + C^". (IV.37) 

where C^" can be determined by comparing r^'^^jj from the effective field theory (|IV.12p with the total energy- 
momentum tensor from the hydrodynamic side (jlV.lOp . It turns out to be arbitrary. To show this, we will 
define, as in [12], an energy- momentum tensor of the fiiud TiCum. This will done for the effective field theory, as 
well as for hydrodynamics. Only under the assumption that these two tensors are equal, it can be shown that 
C^"' in PV.37P vanishes. 

Let us start by defining, as in , the energy-momentum tensor of the fiuid Tfi^]^ ^'^ using the notation of 
the hydrodynamics 

T^^,/" = wu'^u'' - Pn^.^g'"' - BMu^'u", (IV.38) 

where Pn„M = -V^^^ = P + V^^\ with P = -17 = + V'^^^). Here, and V^^^ are the tree level and 

one-loop effective potentials, respectively. They are defined in pi.7[) generally and for our specific model in 
(IIL221) . In (irOSll . rCd satisfies 

u^d^T^d,/"" = u^Rj^d^'pk - u^d^'Vo - u^a^^i, (IV.39) 

where pV.34p and the definition Rk from pil.3[) are used. On the other hand, let us consider the energy- 
momentum tensor pV.12p of the effective field theory, and define^^ 

ri')/" = Tk + ^l.g^^ (IV.40) 

where Tk is defined in pV.13p . Using the EoM, we have already shown that df^Tj^" = 0. We arrive therefore at 

u^d^T^^^l^" =u^d^'Vi. (IV.41) 

Comparing pV.4ip with pV.39|) and assuming that T^l}^^" — Tf['^}^'^'^ , we arrive at the relation 

Tk = {e + P-BAI)u''u'', (IV.42) 

where Pf,„M ee ~V^^^ = P + and P = -Q = -(y(o) + V^^^) are used. Plugging finally pV.42p back in 
PV.12[) and comparing the resulting expression with pV.37[) . we get C^" — 0, and 

^(tot) = i^ + P- BM)u^u'^ - F^^M^\ with u,d,^:,^ = u^KJitot), 

as expected. Note that the above assumption Tf^^^'^" — Tal}^^" leading in particular to the relation pV.42p 
is the link between effective field theory and hydrodynamics. Whereas Tk, on the effective field theory side 
of pV.42p . includes the detailed information concerning the matter content of the model and the interactions 
involved, the expression {e + P — BM)u^u'' on the hydrodynamic side describes effectively the state of a perfect 



See [14| for a similar definition. 
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magnetized fluid through its energy density e, its pressure P and its magnetization M. In contrast to [4J|, 
the present MHD description of a magnetized fluid consists of chiral field or equivalently quasiparticles a and 
vr in the language of condensed matter physics, and exhibits through the special construction of our model a 
spontaneous chiral symmetry breaking. 

In the next section, using the above MHD description and performing a first order stability analysis, we will 
study the effect of the external magnetic field on the sound modes propagating in an expanding magnetized 
QGP including the chiral fields. 



V. STABILITY ANALYSIS AND SOUND VELOCITY IN A STRONG MAGNETIC FIELD 



In this section, we will first determine the dispersion relation of a perfect magnetized QGP at finite (T, /i) 
coupled to chiral fields. Using the dispersion relation, we will then determine the sound velocity of plane waves 
propagating in this hot and dense medium. Due to the broken rotational symmetry in the presence of a constant 
magnetic field, we expect the sound velocity to have an unisotropic distribution. In particular, for plane waves 
propagating in the transverse as well as longitudinal directions with respect to the external magnetic field, 
different dependence on the angle (0, Lp) of the spherical coordinate system will arise. 



A. The dispersion relation in a perfect magnetized QGP coupled to chiral fields 



To study the effect of the magnetic field on the magnetized plasma modeled in this paper via the effective NJL 
model in the presence of a strong magnetic field, we will study the onset of instabilities by performing a first 
order stability analysis. In Table [H a list of all the necessary relations from field theory, thermodynamics and 
hydrodynamics is presented. They are assumed to be valid in a state out of equilibrium. 



Field theory 




= -Rk, from (|III.2|) 


Thermodynamics 


e + P = Ts + p,n + BM, 


from PV.6|) 




de = Tds + pLdn + BdM + Rkdpk, 


from PV.12|) 




dP = sdT + ndp + MdB - Rkdpk, 


from (1IV.15I) 


Hydrodynamics 


d^inu") = 0, 


from PV.23P 




d^isu'') = 0, 


from (|IV.24|) 




u,^^, {{e + P- BM)u^'u") = u^'d^^P + u^'Rkd^^pk, 


from PV.34|) 



TABLE I: 



As was indicated in the previous section, we will assume that B, T, /i, and M remain constant (in the space-time 
x) and will expand n, s, p = (cr, tt), and around their equilibrium configurations no, sq, po and Vq = (1, 0) to 
first order 

n{x) = no + ni{x), 
s{x) = Sq + si{x), 
a{x) = (Tq + (Tiix), 

Tr{x) = TTo + TTl{x), 

u^ix) = (V.l) 

Here, ttq = 0, Vq — (1,0) and = (0,vi). Further, according to the arguments in the paragraph following 
(jIV.lTp . in the thermal equilibrium R^g = i?7r„ = 0. Perturbing first the EoM (jIII.2[) . we arrive at 

Pra.aVi,. + ^^Pr^.aVi. = -ni (^) (^) - ml PI,, (V.2) 
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with the mass matrix 

(V.3) 



'■ke 



Here, we have used the definition of Rk = which we assume to be also valid in a state out of equilibrium. 
Expanding e(p; n, s, M) up to second order in the thermodynamic quantities, we get 

with Co — eo {p'o^no, sq, M). Differentiating (|V.4p with respect to po,fc and using i?cro = Rno — 0; we arrive at 
fyT2)) . Note, that Fl'{ao;eB,T, fi),i = 1,2 are already defined at the vacuum configuration (ctq 7^ 0, ttq = 0) 
[see Sect. HTC] . For fc = tr, (fO]) reduces to 

[G^'a.a- + .4] = -n, f^) - (^) , (V.5) 

with TO^^ = and G''^(o'o; eS, T, fi) = + 2F^, whereas for fc = tt, we arrive at 

FCa^a^TTi = 0. (V.6) 
The conservation laws (|IV.23p and pV.24p lead to 

9oni(a;) + rioV • vi — 0, 

5oSi(a:) + soV-vi = 0. (V.7) 

Finally, perturbing the energy- momentum conservation law pV.34[) . we get 

W^oV-vi = -^doei, (V.8) 
Wo 9ovi = -VPi, (V.9) 

where Wq = eo + Po - BM ^Ts + pn. Using 

ei=^ni+Tsi, (V.IO) 

that arises from (jV.4[) . and the conservation laws (jV.7p . the time component of the Euler's equation, (jV.Sp . is 
automatically satisfied. As for ()V.9p . Pi is given by 

dPo ^ dPo ^ dPo 

-Pi="-i^ hsi^ hcTi-^ — , (V.ll) 

ono oso oao 

where Pq = ^'(eoj "-Oj so, ctqi -^)- Using now Pq = — eo + Tso + fJ-no + BM, and the thermodynamic relations 
(jV.14[) in thermal equilibrium, wc have 

no^r^ + So- — 7^ h M- 



dno driQ dsodno dMdno 

dPo d\o , 52 £0^,, d^eo 

= no- — h so^^TT + 



9(70 
with 

1 



dso dnodso Osq dMdso 

dP 

° K [Wo + BM) = R'^ (eo + Po) , (V.12) 



[Wo + BM) 



no 



(V.13) 
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To determine the dispersion law in the above magnetized QGP, we consider a plane wave of the form ^i(x) — 
^^^-ikx^ with ^1 = {ni, si, (Ti, TTi, vi} and fc^ = (w, k) in the medium and rewrite the relations (|V.5P - (jV.7p . and 
(|V.9P in the momentum space as 



0, ,2 



.2 1,2 



ml) CTi 



0, 
Wo 



(tj2 - k^P') k • Vl = 



— iWo + BM)k'Kau 

Wo 



(V.14) 



where, 



K 
p' 



1 






Wo 


no 


y dn 


1 


no 


rdp\ 


Wo 





So 



So,M,tTo 



so,M,ao 



So 



\ ds 

dP 
97 



no,M,ao 



no,M,ao 



In Appendix [Cl we present a derivation of the last two equations in (|V.14[) . using the relations 

k • vi 



no 



Sl 

So 



(V.15) 



(V.16) 



that arise by plugging ni{x) — fiie as well as si(x) 
in (|V.14p . we arrive at the dispersion relation 



sie '^''^ in (|V.7p . Multiplying the last two equations 



ml) (lo^ 



- k^P') = k^R^Ri (Wo + BM) . 



(V.17) 



In what follows, we will linearize the above dispersion relation in k and determine the sound velocity of a 
plane wave propagating in the magnetized QGP coupled to the chiral field a. Note that, in the above first 
order perturbation, the pion field tt is decoupled from the dynamics of the expanding magnetized QGP. Similar 
phenomenon is also observed in [l^ , where the effective action of the linear a model is used to present a chiral 
hydrodynamic description of an expanding QGP. 



B. Sound velocity in a perfect magnetized QGP coupled to chiral fields 

Sound is defined as a small disturbance propagating in a uniform, field free fluid at rest. Perturbing the energy 
density e(x) and the pressure P{x) around their equilibrium values eo and Pqj up to small Se and SP, and 
linearizing the energy-momentum conservation equation, we arrive at 



de 



ciAiSe) = 0, (V.18) 



where = For an ideal classical gas, consisting of independent and massless particles without interaction, 
Cg = 1/3. This can be shown using the identity P = |, which is derived using the kinetic pressure and the energy 
density in terms of Maxwell-Boltzmann statistics, that replaces the Bose-Einstein and Fermi-Dirac statistics in 
a classical limit. Note that the relation P = I, leading to — 1/3, arises only by assuming the rotational 
symmetry. It holds approximately for a uniform QGP at high temperature, where interaction arc small due to 
asymptotic freedom |l|. 

In [l3 | , the sound velocity of a perfect chiral fluid is determined using the effective potential of a linear a 
model in terms of chiral fields. The dispersion relation of this model is similar to (|V.17[) . where G* = 1, « = 0, • • • 3 
and B = 0. Assuming the rotational symmetry, the dispersion relation is then solved. For long wavelength 
fluctuations, the roots are approximately given by (l^ 



2 2 

uj„ = m„ 



0{k'), 

^l/k' = (p'-^]+0(k'), (V.19) 



23 



where Wq = eo + Pq. Note, however, that m the presence of a constant magnetic field which is ahgned e.g. in the 
third direction, the rotational symmetry is a priori broken, and the sound velocity will be therefore different 
from = 1/3 of an ideal non-interacting QGP. It will be also different from (|V.19p for a perfect QGP coupled 
to chiral fields. 

To determine the sound velocity in a perfect fluid in the presence of a constant magnetic field and consisting 
of chiral fields, we will use a method introduced in (l6j |. where a linearized hydrodynamical theory of magnetic 
fluids^^ in a strong magnetic field is presented. Solving the equations for a sound wave propagating at angle 
{9, (f) with respect to the external magnetic field direction, the sound velocity is shown to be anisotropic. Here, 
{9, if) are the angle in spherical coordinate system. In the next paragraphs, we will consider two different cases 
separately: a) propagation in the ei — 62 plane, and b) propagation in the ei — 63 plane. We will show that 
whereas the the sound velocity for a propagation transverse to the external magnetic field B = Be^ (case a) is 
independent of the angle ip, for the propagation longitudinal to B (case b), it depends on the angle between 
the wave vector k and the direction of the external magnetic field. 

1. Propagation in ei — 62 plane 

Let us begin our derivation, by choosing a plane wave with the wave vector k = (fe sin 9 cos ip, k sin 9 sin (p^ k cos 9) 
in the ei — 62 plane, transverse to the direction of the magnetic field B = Be^. This corresponds to 9 = 7r/2 
and iy9 ^ 0, where < 9 < -k and < 1^9 < 27r are angles in spherical coordinate system. The wave vector 
is therefore given by ki2 = (fci2 sint/j, fci2 cos(/?, 0), where the superscripts correspond to a propagation in the 
ei — 62 plane. For a given frequency equation (jV.17|) determines k. In a linearized hydrodynamics, we will 
evaluate k around fco = ^ , where Cs = 1/ -\/3 is the sound velocity in an ideal classical gas. Thus plugging, the 
relation 

k^2^ka + k[^, (V.20) 

in the dispersion relation (jV.17[) and expanding the resulting expression in the order of fcjj to 0{k'i2), we 
get 

= ^, (V.21) 

with the numerator 

M2 = -2cl[rnl{cl~P') + R„Kw,]+u:\2clG^-G^-G'){cl-P') 

-cj2(g1-G2)(c2-P')cos(2(^), (V.22) 

with Wo = Wq + BM = eo + ^0 and the denominator 

2?i2 = 2c2[m2(c2 + P')-<-R>o]-w'[2c^G°-3(Gi + G2)P' 

+cl (G^ +G^ + 2G"P')] - (G^ - G^)u^ cos{2ip). (V.23) 

Setting G' = l,i = 0,1,2, P' — and P = as well as R'^ = i?" = 0, we get k[2 ~ 0, as expected. To 
determine the frequency dependent anisotropy A, we use the identity u ~ fci2i'i2 with fci2 — ^ + k'12 and 
Vi2=Cs{l + A12). Using k[2 from (|V.2ip - (|V.23p . we arrive at 

Ai2(G°, G\G^; ^, = -TT^i^, (V.24) 

K12CS + ^ 

where G*,i — 0,1,2 are fimctions of eB,T and fi. In our specific model with G^ — G^ [see Sect. Ill C) . the 
(/7-dependence in the anisotropy A12 vanishes. The anisotropy is then given by 

. .^0^1 . ^ -2c^K(c^ - P') + + c.^(2egG" - 2G'){cl - P') 

i2[u ,u , ) 2c2[m2(c2+P')-fi;r-R>o]-c^2[2c4GO-6GiF' + c2(2Gi + 2GOF')]' 
Setting G* = 1, « = 0, 1, P' = cj. and S = as well as R'^ ~ P" — 0, we get A12 — 0, as expected. 



In [16| . a magnetic fluid is defined as a coUoid of tiny (100 A) magnetic particles or grains suspended in a carrier fluid as water. 
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2. Propagation in ei — 63 plane 



Setting ^ 0, = 0, the wave vector in the ei — 63 is given by ki3 — (fci3 sin 6*, 0, fcia cos 6'). Plugging now 
the approximation ki-^ = fco + k'l^^, with /co = ^ in the dispersion relation (|V.17p . and expanding the resulting 
expression in the orders of fc5^3, we arrive at 

fc'i3 = (V.26) 

with the numerator 

M3 = 2c^c.K(c2-P') + <i?>o]-^'(2c^G°-Gi-G3)(c?-P') 

-w3(Gi - G3)(c2 - P') cos(20), (V.27) 

and the denominator 

P13 = 4cV(Gi+G3)P' + 2c^[2m2p'_2i?^i?>o-w2(Gi+G3 + 2G"P')] 

+2c,w2(Gi - G3)(c2 - 2P') cos(26i). (V.28) 

Setting G* = l,z = 0, 1,2, P' = and P = as well as P^ = P" = 0, we get ^'13 = 0, as expected. The 
anisotropy A13 is then determined using the identity lj « fci3Wi3 with fci3 = fco + ^13 and t;i3 = (1 + A13). 
Plug ging k'lQ^ from (jV.26[) - (jV.28[l in fci3, the anisotropy reads 

Ai3(GO,Gi,G3;0,t^) = ^, (V.29) 



with 



-2cl[ivl(cl P') + p:p>o] + u\2clG^ - - G'){cl - P') 

+oj^{G^ -G^){cl~P')cos{2e), (V.30) 



and 



Ad = 2cl[ml{cl + P') - P;,P>o] - c^'[2c^G0 - ?,{G^ + G3)P' + c2(Gi + G^ + 2G°P')] 

W(Gi - G3)(c2 - 3P') cos(20). (V.31) 

Setting G^ ^ \.i ~ 0,1, P' = and P = as well as P^ = P" = 0, we get A13 — 0, as expected. Here, 
in contrast to A12, the ^-dependence remains in the anisotropy A13. In the next section, we will determine 
numerically the anisotropy functions A12 and A13 for plane waves propagating in the transverse and longitudinal 
directions with respect to the external magnetic fields, respectively. 



VI. NUMERICAL RESULTS AND DISCUSSION 



A. Energy density and pressure of a magnetized QGP near tlie chiral phase transition point 

In this section, we will use the results from previous sections to study the effect of a strong magnetic field on 
the sound velocity Vs of a plane wave propagating in an expanding magnetized QGP. The latter is modeled in 
this paper by the NJL effective action in the presence of a constant magnetic field. In previous section, using 
a linear approximation, we have defined an anisotropy function A = ^ — 1, where Cg = l/-\/3 is the sound 
velocity in an ideal gas without matter fields. In js^l, the temperature dependence of Vs in a hot QGP without 
magnetic field is studied. It is shown that the sound velocity shows a minimum at the chiral critical point. 
Whereas the results in [s^] is found using purely hydrodynamical arguments, we intend to present in this paper 
a microscopic model which is merged carefully with hydrodynamics to study the (T, /i) dependence of sound 
velocity in a magnetized QGP. To this purpose, we have to know approximately at which critical temperature, 
Tc, the expected chiral phase transition may occur. The temperature dependence of the energy density of a 
system can give us a useful hint in this direction. As it is known from lattice QCD, we expect an increase in 
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pressure P and the energy density e at the chiral critical point as a consequence of a first order phase transition. 
To determine P in our NJL model in a strong magnetic field, the relation (jV.lOl) can be used between the 
thermodynamic pressure and the one-loop effective potential which is determined in Sect. Ill Al The energy 
density e can be determined using the equation of state (jV.6[) . i.e. e = —P + nfi + TS + BAI. Here, the number 
and entropy densities n and s, as well as the magnetization density M are given by thermodynamic relations 
(|V.16p . To determine the thermodynamical quantities e,P,n and s at thermal equilibrium, i.e. to,PQ,nQ and 
So, we have to replace p — (o", tt) appearing in the thermodynamical potential (|II.22p by the dynamical mass 
ao{eB,T, ii) which is determined in (jll.29|) .^^ As we have shown in Sect. IIIBl the dynamical mass can only 
be determined approximately. Using the same approximations as described in the last paragraph of Sects. Ill Al 
and III B1 we have determined the energy density eq/T'^ and pressure 3Po/T^ as a function of T for a fixed value 
of the magnetic field eB = 1 GeV^ and two different values of the chemical potential fJ, — 0, 10^^ GeV in Fig. 
21 As it is expected, the energy density and pressure have a maximum at Tj, ^ 100 — 120 MeV. 



a) eB = 1 GeV- and // = GeV 



b) eB = 1 GeV^ and fi = 10"' GeV 




0.1 0.2 0.3 0.4 0.5 0.6 
Tin GeV 




- 3 Po/T" 



0.1 0.2 0.3 0.4 0.5 0.6 
T in GeV 



FIG. 4: The temperature dependence of the energy density eo and pressure Po for eB = 1 GeV'^ and /i = GeV (a) as 
well as ^ = 10'^ GeV (b) at thermal equilibrium. At Tc ~ 100 - 120 MeV, there is an increase of the energy density as 
a consequence of a first order phase transition. 

We will next study the effect of the strong magnetic field fixed at eB — 1 GeV^ on the sound velocity in the 
vicinity of the chiral critical point Tc ~ 100 — 120 MeV. As we have seen in Sect. IV B ll and lVB2| the presence 
of an external magnetic field leads to an anisotropy in the sound velocity, so that the anisotropy function A for 
a plane wave propagating in the transverse (A12) and longitudinal direction (A13) with respect to the external 
magnetic fields have different dependence on the angles {9, ip) of the spherical coordinate systems. In what 
follows, we will study the temperature dependence of A12 and A13 separately. 



B. Anisotropy function and sound velocity for a propagation in ei — 62 plane 

As we have seen in Sect. IV B II for a plane wave propagating in ei — 62 plane, transverse to the magnetic field 
B = Be^, the anisotropy A12 from (|V.25|) does not depend on the angle ip of the spherical coordinate system. 
On the other hand, A12 depends on the frequency of the incident plane wave oj. Fig. [5] shows A12 as a function 
of OJ for fixed values of eB = 1 GeV^ and ^ = GeV and for T = 50 MeV (below T^),T = 100 MeV (- 
and T = 150 MeV (above Tc). The three curves show different slopes, but the maximum change of A12 for the 
frequency lo e [0, 10] fm~^ is small ~ 0.2%.-*^^ In what follows, we will fix the frequency io he uj = 0.1 fm~^, 
where according to Fig. [5l A12 remains approximately constant for T > Tc. 

In Table [III the anisotropy function A12 and the velocity V12 are determined for fixed eB = 1 GeV^ and various 
T S [30,450] MeV and /i = 0, 10^'^ GeV. Whereas for /i = 0, the anisotropy A12 and the velocity V12 are real 
valued, they are imaginary for /i 7^ 0, e.g. /i = 10""^ GeV. Qualitatively, the real part of A12 as well as V12 for 



Let us remind that the configuration (o-q ^ 0, ttq = 0) builds the stationary point of the effective potential or equivalently 
describes the thermodynamic equilibrium. 

In [38| | for the relevant frequency is taken in the interval u) 6 [—0.1, +0.1] fm~^. 
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FIG. 5: The frequency dependence of A12 = ^ - 1 for three different temperature T = 50 MeV (below T^), T = 100 
MeV (~ Tc) and T = 150 MeV (above T^) for = 1 GeV~^ and zero chemical potential. 



^ = and /X 7^ has a maximum at T 40 - 45 MeV~ 0.4 - 0.45Tc. In the transition regime ~ 100 - 150 
MeV, they arrive at their local minimum and remain constant after phase transition T > 150 MeV (see Fig. |6]). 
The constant values of Vyi after the phase transition is = 1-5 Cg 0.866. The temperature dependence of 
the sound velocity is studied in 37| for eB — and /i = using purely hydrodynamical methods: The authors 
use at high temperature T > 1.15Tt, the sound velocity function obtained from lattice QCD [s^, whereas at 
low temperature T < 0.15Tc the result of the hadron gas model is used. At moderate temperatures different 
interpolations between those two results are employed. It is shown that the interpolating functions have a local 
maximum at T = 0.4rc (corresponding to the maximal value of the sound velocity in the hadron gas) and a 
local minimum at T = (corresponding to the expected minimal value of the sound velocity at the phase 
transition). Comparing to the results from 37|, our results show that for a wave propagating in the ei — 62 
plane, transverse to the direction of the magnetic field, the qualitative dependence of the sound velocity v\2 on 
temperature does not change. As for the imaginary part of vyi for ^ 0, as it can be seen in Table HIl they are 
several orders of magnitude smaller than the real part of vn- These kind of density fluctuations are studied in 
38[, and are interpreted as a possible origin for the suppression of Mach cone at the chiral critical point. 



eB = 1 GeV^ ^ = GeV, uj = 0.1 



eB = 1 GeV% M = lO"-^ GeV, u) = 0.1 



T in MeV 



Ai 



f 12 



Ai 



1^12 



30 
40 
45 
50 
100 
150 
200 
250 
300 
350 
400 
450 



0.518087 
0.54583 
0.569558 
0.513846 
0.500001 
0.500001 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 



0.876468 
0.892486 
0.906184 
0.874019 
0.866032 
0.866026 
0.866025 
0.866025 
0.866025 
0.866025 
0.866025 
0.866025 



0.516894 + 2 
0.538411 - 6 
0.522862 + 1 
0.528016 - 1 
0.500005 - 3 
0.5-2.67 
0.5 - 4.98 
0.5 - 2.46 
0.5-2.97 
0.5-8.57 
0.5 + 8.75 
0.5-1.52 



.87 X 10"^ 
.81 X 10"^ 
.17 X 10"^ 
.69 X 10"^ 
.31 X 10"^ 
X 10"^**i 
X 10-3*i 
X 10~^^i 
X 10~^"i 
X lO^^^i 
X lO'^'^i 
X 10-^^ 



0.875779 
0.888202 
0.882020 
0.882200 
0.866029 
0.866026 
0.866025 
0.866025 
0.866025 
0.866025 
0.866025 
0.866025 



- 1.66 X 

- 3.93 X 

- 3.93 X 

- 9.73 X 

- 1.91 X 

- 1.54 X 

- 2.86 X 

- 1.42 X 

- 1.72 X 

- 4.95 X 
+ 5.05 X 

- 8.75 X 



^-28 
-31 
1-31 
33 
■27 



10 

^Q-28 



-34 
-33 

10-3° 
10-31 

10-35 

— ■51, . 



10" 
10' 



10-32 



TABLE II: The anisotropy function A12 and the velocity W12 for T G [30, 450] MeV and fixed eS = 1 GeV^ For nonzero 
baryon chemical potential A12 as well as 1112 are imaginary. At the transition temperature A12 as well as v\2 reach their 
minimum value and remain constant after the phase transition T > Tc ~ 100 — 150 MeV. 
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a) eB = 1 GeV^ /i = GeV and m = 0.1 fm"' 



b) eB = 1 GeV^ /i = 10"' GeV and oi = 0.1 fm"' 
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FIG. 6: The anisotropy A12 = ^ — 1 as a function of temperature T for fixed value of eB = 1 GeV^, = 0.1 fm ^ 
as well as /I = GeV (a) and fi = 10"'' GeV (b). The sound velocity has a maximum at T ~ OATc — 0.45rc and 
decreases at the critical temperature Tc ^ 100 — 120 MeV. After the phase transition at T > 150 MeV, it remains 
constant V12 ~ 1.5cs. 



C. Anisotropy function and sound velocity for a propagation in ei — 63 plane 

For a wave propagating in the ei — 63 plane, parallel to the external magnetic field, the anisotropy function A13 
from (|V.29|1 as well as the velocity wia = Cs(l + A13), depend on the angle 6 of the spherical coordinate system. 
It depends also on the frequency u>, as for A12. Fig. [7] shows the {u!,6) dependence of A13 for fixed values 
of eS = 1 GeV2, ^ = as well as T = 50 MeV (below Tc), T = 100 MeV (- T^) and T = 150 MeV (above 
Tc). Qualitatively, whereas A12 increases with oj (see Fig. [5]), Ai3(6' 7^ ^) decreases with u; and depends at the 
same time on the angle 9 between the wave vector k and the external magnetic field B. Note that primarily for 
9 ~ ^, the anisotropy is to be calculated from A12 (jV.25|) . But, it turns out that 

Ai2(G°, G';u;)^ Ai3(G°, G\ G^; ^ = |, c.), 

where A12 is given in (jV.25p and A13 in (jV.29|) . Thus, A13 (0 = ^) has the same w-dependence as A12. 



a) Ai3 for fixed eB = 1 GeV^ T=50 GeV and //=0 b) Ajj for fixed eB = 1 GeV^ T=100 GeV and ii=a c) Au for fixed eB = 1 GeV^ T=150 GeV and ji=0 




FIG. 7; Qualitative dependence of A13 on the frequency lu and the angle 6 for fixed values of eB — 1 GeV^, /i = as 
well as (a) T = 50 MeV (below Tc), (b) T = 100 MeV(~ Tc), and (c) T = 150 MeV (above Tc). In contrast to A12, A13 
as well as the sound velocity for ^ 7^ decreases with uj. For 9 ~ A13 has the same behavior as A12. 



In what follows, we will work with a fixed value of = 0.1 fm^^, where for T > Tc, A13 is approximately 
constant (see Fig. [7]). In Fig. [8] a) the temperature dependence of the sound velocity t;i3 is presented for fixed 
values of 6* S {0, f , f , f } and eB = 1 GeV^ as well as /j, = GeV and w = 0.1 fm^^. In contrast to V12, the 
sound velocity W13 increases with temperature. It has therefore a local minimum at T < Tc, reaches its local 
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maximum at T ^ T^, and remains constant for T > T^. Hence, the constant value W13 = 1.5 Cg « 0.866 of the 
sound velocity that arises in the above linear approximation, is a lower bound for V12 and an upper bound for 
W13 as the temperatures are higher than the critical temperature T > T^. For U13 this behavior is independent 
on the angle [see Fig. [5]b for a comparison between V12 and U13 at various angles]. 



a) vi3 for eB = 1 GeV^ // = GeV and o) = 0.1 fm"' 



b) V12 versus V13 for eB = 1 GeV^, /i = GeV and cj = 0.1 fm ' 
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FIG. 8: a) Temperature dependence of the sound velocity V13 for fixed values of 6 G {0, and eB — 1 GeV^ as 

well as = GeV and co = 0.1 fm~^. b) Temperature dependence of «i2 versus V13 for 6' G {f, f , f} and e_B = 1 GeV^ 
as well as /i = GeV and lo = 0.1 fm~^. The sound velocity V12 (via) has a local maximum (minimum) at T < Tc, 
reaches its local minimum (maximum) at T ~ Tc and remains constant for T > Tc- 
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FIG. 9: The ^-dependence of 7J13 for T = 30, 45, 50 MeV (< Tc) T = 100, 150 (> Tc) and for fixed values of eS = 1 



GeV^ /i = GeV and to = 0.1 im'^ Whereas for T < Tc. 
T > Tc it is almost constant with V13 — 1.5 Cs ~ 0.866. 



«13 is 6 dependent and has its maximal value at 9 



for 



In Fig. M the 9 dependence of ^13 is plotted for T = 30, 45, 50 MeV (< T^) T = 100, 150 (> T^) and for fixed 
values of eB — 1 GeV^, /i = GeV and oj ~ 0.1 fm^^. As expected from Fig. [3 whereas for T < V13 is 9 
dependent and has its maximal value at 6 = ^, for T > Tc it is almost constant with vi^, — 0.5cs w 0.866. 

VII. CONCLUSION 



Strong magnetic fields play an important role in the physics of non-central heavy ion collisions. The y p rovide a 
possible signature of the presence of CP-odd domains in the presumably formed QGP phase 23, 2^ 2J|- Apart 
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from other mechanisms, relativistic shock waves, that are beheved to be built in the heavy ion coUisions, can 
be made responsible for their generation. In astrophysics, the strong magnetic fields generated in collisionless 
relativistic shocks play an important role in the fire-ball model for Gamma-ray Bursts [29]. Recent studies on 
the effect of magnetic field in modifying the nature of the QCD chiral phase transition in a linear a-model 
indicate that for high enough magnetic fields, comparable to the ones expected to be created in the non-central 
heavy ion collisions at RHIC, the original crossover is turned into a first order phase transition 3^, 34, 3^. 

In the present paper, we have studied the possible effects of strong magnetic fields on the hydrodynamical 
signals of an expanding perfect magnetized QGP coupled to chiral fields. In particular, performing a first order 
stability analysis in a chiral magnetohydrodynamic framework, the sound velocity of a propagating plane wave 
in this medium is determined. To this purpose, we have extended the variational method applied in 1^ on the 
effective potential of a linear cr-model at finite (T, /i) and zero magnetic field B to the case of nonzero magnetic 
field. As for the effective action, we have used the NJL model of QCD at finite (T, /i), and in the presence of a 
strong B field. The magnetic field is assumed to be constant and aligned in the third direction, i.e. B = Be^. 
This model is known to exhibit a dynamical chiral phase transition due to the phenomenon of magnetic catalysis 



arising from the dynamically generated fermion mass 19 1 



In the first part of this paper, we have explicitly determined the effective potential and effective kinetic term 
of the NJL model at finite (T, /i) and in the presence of a strong B field. In particular, the effective kinetic term 
is determined in a derivative expansion, where structure functions appear that depend on (T, /z) and the external 
magnetic field B. Using the effective Lagrangian density, we have derived the dynamical mass generated in the 
regime of LLL dominance. The dynamical mass, being the configuration that minimizes the effective potential 
of the model, is used, in the second part of the paper, as the equilibrium configuration once the instabilities in 
the strongly magnetized QGP are set on. It depends on (T, /i) and B and vanishes once the QCD matter has 
passed the chiral critical point (see Fig. 2 for its qualitative behavior for a fixed value of the magnetic field). 

In the second part of the paper, a chiral magnetohydrodynamic description of the strongly magnetized QGP is 
presented. We have, in particular, compared the total energy- momentum T^"^ and the polarization tensor Af'', 
that arise field theoretically from the Lagrangian density of the effective NJL model, with the corresponding 
quantities from hydrodynamics. The polarization tensor M^'' appears in T^" and includes the magnetization 



density M of the medium. It is introduced in [4J], where Nernst effect near the superfluid-insulator phase 



transition of condensed matter phy sics is studied using an appropriate MHD description. Note that in contrast 



to the work by Sachdev et al. |44| . the magnetized fluid modeled in our paper is coupled to chiral fields and 



involves the effects of a possible chiral phase transition. 

The dispersion relation of the hot and dense medium in the presence of a strong magnetic field is derived 
using a first order stability analysis. As for the sound velocity of a plane wave propagating in the magnetized 
QGP modeled in this paper, we have used a method presented in 16| and determined the anisotropy A in a linear 
approximation. Here, A = ^ — I with Vg the sound velocity in the medium, and Cg = I/\/3 the sound velocity 
in an ideal gas. In our model, A shows, as in the magnetic fluid considered in [l^ . a non-trivial frequency 
dependence. Moreover, for a plane wave propagating in the transverse ei plane with respect to the external 
magnetic field B = Be^, A12 does not depend on the angle (p between the wave vector k and the external 
magnetic field B. On the other hand, for a plane wave propagating parallel to the magnetic field in the ei — 63 
plane, A13 depends on the angle between k and B. Here, {0,(p) are the angles in the spherical coordinate 
system. In Figs. [5]b we have compared the sound velocity U12 and U13 as a function of temperature for a fixed 
values of eB, fj,, to and various angles 9. In contrast to the sound velocity W13 of a plane wave propagating in the 
ei — 63 plane, the sound velocity in the ei — 62 plane, V12, has a local maximum at T < Tc, reaches its local 
minimum at T ~ Tc and remains constant for T > T^. The constant value Vg ^ 1.5 Cg « 0.866 seems therefore 
to be, in this linear approximation, a lower bound for V12 and an upper hound for U13 as the temperatures are 
higher than the critical temperature T > Tc- For vi^ this behavior is independent on the angle 9. 

The universal properties of the equation of state e(p), and in particular the temperature dependence of the 
speed of sound at T '--^ (2 — 3)rc, where the theory is close to being conformal, is recently investigated in 54, HH] 
using the holographic principle in the framework of gauge-gravity duality of string theory. It is shown, that for a 
general class of strongly interacting theories at high temperatures the speed of sound approaches the conformal 
value Cg = f/V3 from below. It would be interesting to extend the method used in these papers to study the 
behavior of the sound velocity in a strongly magnetized QGP. 
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APPENDIX A: HIGH TEMPERATURE EXPANSION AND BESSEL FUNCTION IDENTITIES 

FOR FINITE T AND ^ 

As we have seen in Sect, ni the temperature dependent part of the effective potential as well as the gap equation 
of the NJL model of QCD consists of terms in the form 



oo 

5^(-l)^-Xp(fe)cosh(M/?£). 



(A.l) 



Here, z = to/3, with to the dynamical mass,^*^ and /i the finite chemical potential. This expression can be 
written in the generalized form 



C30 

5]-ifp(fe)cos(£0), 



(A.2) 



with = TT + This identity is determined for /i = in 40] . To derive the corresponding Bessel function 
identities for summations of the form of (|A.2p , the authors starts in particular with the identity 



z COS' 



7 + E 



+ (0 - 27r^)2 



27r. 



(A.3) 



where the notation is used to indicate that singular terms in the summation over all values of ^ g (— oo, +oo), 
here £ = in 1/|^|, are omitted ^40]. In this appendix, we will first evaluate the sum over I on the r.h.s. of (jA.Sp 
and present it in terms of a series expansion in z = to/?. Then, using the result for = tt + as it appears in 
the effective potential and gap equation in Sect, [ill we will determine \A.2\ for p — Q and p — I. The method 
presented in this appendix can be used to give and expansion of (jA.2[) in terms of z for all p. This will be in 
particular useful if we are interested in a high temperature expansion of the effective potential, as well as the 
gap equation. 



1. High temperature expansion of (|A.2[) for p = 

To start let us consider (|A.3p and denote the sum over £ on the r.h.s. by (7(0, z, z'), where z = to/3 and z' — ^(3 



+ C((/.,z,z'), 



^Ko{lz) cos( 
1=1 

where in particular 

C(^,z,z') = f E 
Expanding C(7r, z, z') in the orders of z, we arrive first at 

oo oo oo 

C(^,z,z') = E -4,(^,0, z')+ E B,{^,z,z')^Y.^,, 



^z2 + (tt + iz' - 27r^)2 2tt\ 



(A.4) 



(A.5) 



(A.6) 



1=1 



In the case of the NJL model, m should be replaced by (tq. 
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where 



2\1 ^ 2£ + '-fl' 



is the zeroth order term in z, and 

1 1 



Bi{tt,z,z') = < - 



47r2 |i _ 2£+ 167r4 |l _ 2^ + 327r6 |i - 2£ + 

denotes the higher order terms in z. Further, Ci = Let us first consider the sum over I in Ai 



A > 7 I 2 ^ I 2£ + 1 + Hi 



1 



2^+1 



E 



(2£+l) 



Expanding At in the orders of ^ = we arrive at 



^ A(7r,0,z') = EE 



/c=0 £=0 

OO 



(2^+l)2fc+i 



2fc 



(2£ + 1) 

«=0 ^ ' k=l 



where the definition of the Riemann ^-function 
is used. Using now the identity 



y ^— - y 1 = y i^i^ = in2 

£=0 i=\ i=\ 



we get 



(-l)fe (22fe+i _ 1) /z' 



22fe+l 



2fe 



- C(2fc + 1). 



(A.7) 



(A.8) 



(A.9) 



(A.IO) 



(A.ll) 



(A.12) 



To evaluate Bi on the r.h.s. of (|A.6p with given in (jA.Sp . we spht the sum over I e (— oo, +oo) in a sum 
over t G (— OO, 0] and £ e [1, oo). We arrive at 



J2 Bein,z,z') = J2 



(-1)" /z 



42"+i (n!)' 



(^) 



2n 



as weU as 



£s,(7r,z,z')-E 



(-1)" /z 



«=1 



42»+i (n!)' 



1 iz' 
2^2^ 



1 iz' 
2^2^ 



(A.13) 



(A.14) 



Here, (f>''P\z) is p-th derivative of the polygamma-hmction tjj{z) = ^lnr(z) with respect to z. Plugging 
(|XT2l) - (|XT4)) in (fOl) . C(7r, z, z') is given by 



CKz,z') =ln2 + 5:i^i4i^^^^^ 



- C(2fc + 1) 

TT ' 



y 

6l (n!)' 



(^) 



2n 



^(2«) 



1 iz' 
2^2^ 



^(2n) 



1 iz' 
2^2^ 



(A.15) 
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This leads to the first Bcssel-function identity 



oo ^ 

Ko{£z) cos(£(7r + iz')) = - 



7 + In - 



(-ir 



^ 42n+l („n2 
n— 1 V J 



A;=l 



^ * 2 27r 



^ ' 2 27r 



(A.16) 



where (|A.4p is used. 



2. High temperature expansion of (|A.2[) for p = 1 



In what foUows, we wiU determine (|A.2p for p = 1, and cf) = it + z' with z' = by generalizing the method 
used in [40|. Using the relation 



, oo oo 

— V ^Kiiiz) cos(f(/)) = -z V ^o(^2) cos 



dz^ e 



that follows from the recursion relation 



^KAz)^-K,^^{z)--^K,{z), 



we have 



oo / oo 

^ -ifi(fe)cos(£0) = -- dz iz^Kaiez) 
1=1 •' \ 1=1 



cos 



C(0) 



(A.17) 



(A.18) 



(A.19) 



Here, C{(f>) is an unknown function of 4>, that, in contrast to the derivation presented in [40|, is a function of 
z' = In particular, we will set (j) ~ tt + iz' to evaluate Bessel- function identities in the form (|A.ip . To 
determine C{(j)), we multiply (jA.19[) with z and use the behavior (|A.2ip of the Ki,{z) in limit z ^ 0, 



z 



to arrive at 



oo 



e=i 



where the polylogarithm- function Lip(z) is defined, in general, as 

oo ^ 



k=l 



(A.20) 



(A.21) 



(A.22) 



The second Bessel-function identity can therefore be derived from (jA.lQp in combination with (|A.16p and (|A.2ip 



as 



oo _j 

^-i^i(fe)cos 



e=i 



l-27s-21n - 

TT 



^ (-l)"(l + n) /^^^2« 
(r(2 + n))2 

3 



(2n) 



1 iz' 
2^2^ 



(2«) 



1 iz' 
2^2^ 



_z2>p (-l)'=(22'=+i-l) ^z'^ ''^ 



k=l 



22fe+l 



/ \ 2/i: 1 

M C(2fc + l) + ^{Li2(e*("+'"'^) +Li2(e-'("+'"'))}. (A.23) 



Same method can be used to determine (|A.2p for arbitrary p > 1. 
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APPENDIX B: MELLIN TRANSFORMATION AND THE SUMMATION OVER MATSUBARA 

FREQUENCIES AT FINITE T AND /i 



As we have shown in Sect. Ill CI the integrals arising in the expressions for the structure functions G'',yU = 
0, • • • ,4 from (jTL46j , have the general form 

where, fco = {wi — iy,) with tot = ^(2i?+ 1) the Matsubara frequencies, and ji the chemical potential. In 41], a 
similar integral as in (jB.ip . with a — 0,t — Q and ^, — Q potential is evaluated using the Mellin transformation 



-1 /'C+ioo 

fix) = — / x-'M[f;s]ds, (B.2) 



and its inverse transformation 



oc 



M[f;s] = / x'-'f{x)dx. (B.3) 



As it is denoted in |4l| . the above transformation normally exists only in a strip a < 9le[s] < (3, and the inversion 
contour must lie in this strip a < c < (3. In this appendix, we will present a generalization of the results in 41 1 
for arbitrary a, t and nonzero chemical potential fJ- =/= 0. To this purpose, let us start with fi^e ~ il^) on the 
l.h.s. of (|B.1|) . Using (IB.2|) . we arrive first at 



^ f{uji~ifi) ^ —— ^ / ds{ujt~ifi)~' M[f;s 



where the definition of is used. Performing now the sum over £, as 



2 2tt J ' V 2 27r y ' V 2 27r , 
and plugging this expression on the r.h.s. of (|B.4p . the integral / is given by / = + with 

On the other hand, comparing (|B.3[) with (|B.1[) . M[f; s] can be determined. It is given by 

/f rl'^-h 1.2a 2t 

To evaluate this integral, we combine two vectors k and ytoaD = s + d + 2(t + a)-dimensional Euclidean 
vector. Evaluating now the I?-dimensional integral using standard identities from dimensional regularization 



56|, A1[/; s] is given by 



^^^^5+2* dVLd+2a J 



1 r(| + a)r(f + t)r(a-f -f-t-g) 
2(47r)'i/2 r(|)r(a) m2("-t-«)-«-'i 



(B. 
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where dfl,] = 



r(d/2) 



is used. Plugging now this expression in I± from (jB.6p and evaluating the contour integration 



over s using the residue theorem, we arrive at 

1 T{^ + a) f27T\-^"+''+^^*+''^ 



2(47r)rf/2r(a)/3 T (f ) \ P 



fe=0 



C(2ia + k-t-a)-d;l-'^ 



m(3 



2k 



(B.9) 



Here, fc = 0, 1, 2, • • • ,oo labels the pole of the F-function at s = 2(a — t — a + k) — d. Being a power series in 
(m/3) , the result from (|B.9p can be easily used whenever a high temperature expansion of Feynman integrals at 
finite temperature and density is possible. Note that a similar result is also presented in [s^, where the integral 
in (|B.1[) is evaluate for arbitrary a,t and zero chemical potential /i = 0. 



APPENDIX C: THE DERIVATION OF THE EQUATIONS IN (|V.14p 
1. Derivation of the second equation in (|V.14p 



The first equation in (|V.14p is already derived in (|V.6p . To derive the second equation in (|V.14p , let us consider 
(jV.Sp . that is derived from ()V.2p by setting k — a and using the definition (jV.3P 
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we arrive at the second equation in (|V.14p 

a; 

Let us now consider (jCSP and write it as 



d(R„,s,M,a) d(R„,n,M,a) 

n „ ^ S- 



" Wo { din, s, M, a) 



where the notation (f^)^ = g^^'^j is used. Generalizing now the identity g^^^,^ ^^.^ 



d(A,C) 

d{T.,^..BM) 



(C.l) 



(C.2) 



(C.3) 



(C.4) 



(C.5) 



dA dA 



det 1^ from [57 1 



for a case with four variables, and replacing the Jacobian on the r.h.s. of (jC.SP by , 

we get first 



r / 



dn 


dn 


dn 


dn 




-1 


/ 


dR„ 


dRa_ 


dR„ 


dR, 




dR„ 


dR„ 


dR„ 


dR^ 


dT 


df-L 


dB 


da 




dT 


d^i 


dB 


da 




dT 


dfi 


dB 


da 


ds 


as 


ds 


ds 








ds 


ds_ 


ds 


ds 




dn 


dn 


dn 


dn 


dT 


d^i 


dB 


da 






n 


dT 


d^i 


dB 


da 




dT 


d^i 


dB 


da 


dM 


dM 


dM 


dM 






dM 


dM 


dM 


dM 


— S 


dM 


dM 


dM 


dM 


dT 


d^i 


dB 


da 








dT 


d^i 


dB 


da 




dT 


dfj. 


dB 


da 


da- 


dij_ 


da 


da 


) 




[ 


da 


da 


da 


da 




da 


da_ 


da 


da 


dT 


df-L 


dB 


da 




dT 


d^i 


dB 


da 




dT 


dfj. 


dB 


da 



) , (C.6) 



/J 



35 



where the subscript means that the {n, s, a) in the final result of the determinants are to be replaced by 
(riQ, So, (To) from the thermal equilibrium. At this stage the thermodynamic relations from pV.16p can be used 
to get the final expression for R'^ 
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where Pq = Po{T, fi, B,ao), and the Jacobian is defined by 
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2. Derivation of the third equation in (|V.14[) 



To derive the third equation in (|V.14|) . let us start with l|V.9p . where Pi is defined in (jV.lip . 
above, = ^le"*'""'^ with = {ni, si, ui, vi}, we get first 
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as expected. In what follows, we present the final results for P' and R'l from (jC.10|) and (jV.lSp 
Using the same method leading from (jC.Sp to (|C.7|) . P' from (jClOp can be first given as 
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respectively. 

> . (C.13) 



Using now the thermodynamic relations from pV.16|) to replace no, sq, M, ao by the derivative of Po with respect 
to fi,T, B respectively, we get 
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with the Jacobian J defined in (|C.8p . As next, using the definition of i?^ fi-om (|V.13p . 
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and foUowing the same method as described above, we arrive at 
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where the Jacobian J defined in (IC.Sp . 
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